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We establish compelling evidence for the existence of new quasi-one-dimensional descendants of 
the d-wave Bose liquid (DBL), an exotic two-dimensional quantum phase of uncondensed itinerant 
bosons characterized by surfaces of gapless excitations in momentum space [O. I. Motrunich and 
M. P. A. Fisher, Phys. Rev. B 75, 235116 (2007)]. In particular, motivated by a strong-coupling 
analysis of the gauge theory for the DBL, we study a model of hard-core bosons moving on the 
A^-leg square ladder with frustrating four-site ring exchange. Here, we focus on four- and three- 
leg systems where we have identified two novel phases: a compressible gapless Bose metal on the 
four-leg ladder and an incompressible gapless Mott insulator on the three-leg ladder. The former 
is conducting along the ladder and has five gapless modes, one more than the number of legs. 
This represents a significant step forward in establishing the potential stability of the DBL in two 
dimensions. The latter, on the other hand, is a fundamentally quasi-one-dimensional phase that is 
insulating along the ladder but has two gapless modes and incommensurate power law transverse 
density-density correlations. While we have already presented results on this latter phase elsewhere 
[M. S. Block et al, Phys. Rev. Lett. 106, 046402 (2011)], we wih expand upon those resuhs in 
this work. In both cases, we can understand the nature of the phase using slave-particle-inspired 
variational wave functions consisting of a product of two distinct Slater determinants, the properties 
of which compare impressively well to a density matrix renormalization group solution of the model 
Hamiltonian. Stability arguments are made in favor of both quantum phases by accessing the 
universal low-energy physics with a bosonization analysis of the appropriate quasi-lD gauge theory. 
We will briefly discuss the potential relevance of these findings to high-temperature superconductors, 
cold atomic gases, and frustrated quantum magnets. 



PACS numbers: 71.10.Hf, 71.10.Pm, 75.10.Jm 

I. INTRODUCTION 

One of the recent challenges in condensed matter 
physics has been to understand quantum phases char- 
acterized by singular surfaces in momentum space. The 
canonical example of such a phase is the free Fermi gas 
(or more generally a Fermi liquid), where the singu- 
lar surface is simply the Fermi surface. Ironically, de- 
spite our immense theoretical understanding of the free 
Fermi gas^ and the interacting Fermi liquid,^ a con- 
trolled and unbiased numerical demonstration of such 
a phase in an interacting microscopic model of itinerant 
fermions moving in two or more dimensions is still ex- 
tremely difficult. The main roadblocks in this numerical 
pursuit are (1) the infamous sign problem in quantum 
Monte Carlo simulationd^ and (2) the anomalously large 
amount of spatial entanglement present in states with a 
Fermi surface,^^ thereby rendering recently developed 
sign problem-free tensor network state approaches inad- 
equate with current techniques.l^^'^ Although the most 
obvious and familiar, the free Fermi gas and Fermi liq- 
uid are not the only examples of phases with singular 
surfaces in momentum space. 

Since the discovery of the cuprate superconductors, in- 
terest has emerged in novel two-dimensional (2D) quan- 
tum phases that fall outside the paradigm of Fermi liq- 



uid theory, but that still have correlations with singu- 
larities residing on one-dimensional (ID) surfaces in mo- 
mentum space. Perhaps the most prominent example 
are states with a spinon Fermi surface, so-called "spin 
Bose metals" , which have both a long history in the con- 
text of the cuprateJ^ and also a renewed interest in the 
context of the organic materials k-(ET)2Cu2(CN)3 and 
EtMe3Sb[Pd(dmit)2]2-'^^'^Here, due to the presence of 
an emergent gauge field when going beyond mean field, 
a Fermi-liquid-like quasiparticle description is inapplica- 
ble, thus making fully controlled analytics very challeng- 
ing and further promoting the importance of numerics. 
Not surprisingly though, numerics suffers from the same 
difficulties as before: (1) the sign problem is likely a nec- 
essary condition for realistic parent Hamiltonians, and 
(2) the beyond boundary law scaling of spatial entan- 
glement is likely the same as that encountered for free 
fermions.'^ 

We focus here on a closely related quantum phase, 
the "d-wave Bose liquid" (DBL). The DBL is an exotic 
quantum phase of uncondensed itinerant bosons moving 
on the two-dimensional square lattice first considered in 
Ref. [T71 Like states with a spinon Fermi sea, even though 
the microscopic degrees of freedom are bosonic, the DBL 
too has a set of gapless excitations residing on ID surfaces 
in momentum space, i.e., "Bose surfaces," and so access- 
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ing the DBL in 2D using fully controlled techniques, ei- 
ther numerically or analytically, is no more tractable than 
for the spinon Fermi sea. However, we make progress 
studying the DBL in a controlled way by continuing 
to employ the heretofore fruitful philosophj^^ of 
building a picture of such a 2D phase through a sequence 
of controlled quasi-lD ladder studies. In fact, we believe 
that the very presence of singular surfaces in momen- 
tum space actually renders ladder studies more informa- 
tive and allows us to circumvent some of the usual nu- 
merical and analytical difficulties. When placed on the 
TV-leg ladder, a given 2D phase with a continuous set 
of gapless modes residing on a ID surface should enter 
a phase that is a distinctive multi-mode quasi- ID de- 
scendant of the parent 2D phase, with a number of ID 
gapless modes that grows linearly with N (used inter- 
changeably with Ly in this paper: N — Ly).'^'^ In such 
quasi- ID geometries, both numerics and analytics are 
on much stronger footing: in principle, potential parent 
Hamiltonians can be solved numercially with the den- 
sity matrix renormalization group (DMRG),'2ill22! which 
can then be supplemented with variational Monte Carlo 
(VMC) calculations^ using appropriate projected trial 
wave functions to help map out the phase diagram. De- 
spite the fact that the DBL is a strong-coupling phase 
that cannot be characterized perturbatively in terms of 
the original bosons, the slave-particle approach on the 
ladders accesses the phase in a novel way that is tractable 
via a gauge theory that can be treated using conventional 
bosonization techniques. 

Guided by this line of attack, we aim to accomplish two 
main goals in this paper: (1) establish the existence of an 
exotic metallic DBL phase, i.e., a "Bose metal," on the 
four-leg ladder, which we will argue is both a close rela- 
tive of the proposed 2D phase of Ref. ITTIand a nontrivial 
extension of the quasi-lD descendant DBL phase discov- 
ered on the two-leg ladder in Ref. TF; and (2) elaborate 
on our previous work in Ref. 19, which argued for the 
existence of a novel, but fundamentally quasi-lD, gap- 
less Mott insulating phase on the three-leg ladder. We 
emphasize that while we have already presented high- 
lights of the latter three-leg gapless Mott insulator in a 
previous work,!^ the former four-leg gapless Bose metal 
constitutes a new and very exciting result. These two 
seemingly disparate phases are actually close relatives: 
they share a common parton gauge theory description, 
can be modeled within the same class of projected varia- 
tional wave functions, and, finally, manifest themselves in 
the same microscopic model, all of which we will demon- 
strate in this work. 

The model wave function for the DBL is obtained by 
taking a product of two distinct Slater determinants and 
evaluating them at identical coordinates (Gutzwiller pro- 
jection): 

^'6(ri,...,rArJ = *di(ri,...,rArJ*d^(ri,...,rArJ, (1) 

where Nj, is the total number of bosons. In the language 
of operators, this corresponds to writing the (hard-core) 



boson operator as a product of two fermionic partons: 

6t(r)=4(r)4(r), (2) 

where the hard-core boson Hilbert space is recovered by 
requiring that the densities of the two partons, pd^ (r), be 
the same at each lattice site, i.e., pd^ (r) = Pdii^) = p(r), 
where p(r) is the boson density at site r. Within a gauge 
theory framework, the mean-field picture of the phase 
consists of two independent species of fermions hopping 
on the square lattice with anisotropic hopping; the phys- 
ical Hilbert space is then obtained by strongly coupling 
the di and d2 fermions with oppos ite gauge charge to 
an emergent U(l) gauge field.'i^E^ If one takes the di 
(^2) fermion to hop preferentially in the x (y) direction, 
then the corresponding wave function [Eq. ([ij with "^di 
(^'^2) being a filled Fermi sea compressed in the x (y) 
direction] has a characteristic d-wave sign structurepl 
that is, the sign of the wave function goes through a 
sequence, H 1 — , upon taking one particle around an- 
other and hence the label "d-wave Bose liquid." Perhaps 
more importantly, as alluded to above, this projected 
wave function has power law singularities of various mo- 
mentum space correlators, e.g., the boson momentum 
distribution and density-density structure factor, resid- 
ing on ID surfaces in momentum space.^ These surfaces 
are perhaps the most distinguishing feature of the DBL 
and are crucial to our identification of quasi-lD descen- 
dants of the phase on the iV-leg ladder. A strong-coupling 
analysis^^ of the aforementioned gauge theory motivates 
a simple microscopic parent Hamiltonian which can po- 
tentially harbor the DBL phase. This Hamiltonian con- 
sists of usual nearest-neighbor hard-core boson hopping 
supplemented with an explicit frustrated four-site ring- 
exchange interaction (see also Sec. II C): 



H — Hj + Hii, 
H., = -J (^i&r+A + H.c), 



(3) 
(4) 



Hk = x5](fei6r+i6|:+£+g6r+fi + H.C.), (5) 

r 

with J,K > 0. We focus on this model Hamiltonian, 
the so-called "J-K model," extensively in this paper. 
For K — 0, we expect a generic superfluid with the 
bosons condensed at q = and off-diagonal long-range 
order. It is in the regime where the ring-exchange term 
contributes appreciably to the overall Hamiltonian that 
the d-wave Bose liquid is expected to onset. In the 
strong-coupling limit of the lattice gauge theory for 
the DBL, the relative strength of the ring term to the 
hopping term, K/J, increases with increasing hopping 
anisotropy between the di and d2 partons. Thus, we 
expect the ring term to potentially stabilize the DBL 
phase, and, as in Ref. [18] on the two-leg ladder, we do 
in fact remarkably find evidence for such a scenario on 
three- and four-leg ladder systems. 
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To conclude this section, we now discuss several lines 
of motivation behind this work. Aside from being in- 
trinsically interesting in its own right, we believe the 
DBL phase is potentially relevant to modern-day experi- 
ments and other theoretical pursuits in a number of con- 
texts, including high-temperature superconductors, ul- 
tracold atomic gases, and frustrated quantum magnets. 
We start by focusing on the first case, where we believe 
the ideas behind the DBL can be used to describe the 
charge sector of a particular conducting, non-Fermi liq- 
uid phase of itinerant electrons, i.e., a "strange metal," 
which may possibly be related to the infamous strange 
metal phase of the cuprates. Specifically, we have in 
mind the following scenario. In the spirit of the slave- 
boson treatment of the t-J model,'!^' one can (excluding 
site double occupancy) decompose the electron creation 
operator as a product of a slave boson ("chargon") and 
a fermionic spinon: c\{r) — 6^(r)/j^(r). This leads to a 
gauge theory formulation in which the spinous and slave 
bosons are coupled to an emergent gauge field. While it 
is natural in this context for the spinous to form a Fermi 
sea (see above), the behavior of the slave boson holds the 
key in determining the properties of the resulting elec- 
tronic phase: if the slave bosons condense, a traditional 
Fermi liquid phase is obtained, whereas if they do not, 
we may say that the phase is a "non-Fermi liquid." An 
example of how the latter case can be achieved involves 
further decomposing the slave boson into two fermionic 
partons just as we did for the real boson in Eq. ([2]): 
6t(r) = 4(r)4(r)> where again di (^2) is taken to fill 
a Fermi see compressed in the x (y) direction, and now 
there are two emergent gauge fields needed to enforce 
the physical electronic Hilbert space. This setup sets the 
stage for a theory of a "d-wave metal" phase. While we 
will present work on this phase elsewhere, both in 2D and 
on the two-leg ladder, we believe scaling up the picture 
of the charge sector to many legs (as we do in this pa- 
per) will be an informative endeavor toward an eventual 
theoretical understanding of the 2D "c?-wave metal." 

A direct experimental realization of the DBL can po- 
tentially be achieved in systems of ultracold quantum 
gases. Although it has been proposed that the boson 
ring-exchange model [see Eg. j3|] can be engineered di- 
rectly in a cold atom system,'2Srperhaps the most feasible 
scenario involves engineering a pair of mismatched Fermi 
surfaces in a two-component Fermi gas by introducing 
a hopping anisotropy between the two species.'^Sl Such 
anisotropy can be achieved with a spin-dependent (more 
precisely "hyperfine-state-dependent") optical lattice, a 
setup that w as exp erimentally first demonstrated sev- 
eral years ago.'^SSSl In the noninteracting limit, we then 
have a situation very similar to the mean-field descrip- 
tion of the DBL phase as discussed above, i.e., a system 
of two independent species of fermions characterized by a 
hopping anisotropy. However, a fundamental difference 
is that we are now talking about real fermionic atoms, 
as opposed to fermionic partons modeling real hard-core 
bosons [see Eq. (pi)]. What are the potential phases that 



such a system could enter upon adding local attractive in- 
teractions? This question was first asked by Feiguin and 
Fisher in Ref. dHl and although there are several possibil- 
ities within a BCS treatment, the most interesting sce- 
nario involves formed Cooper pairs entering a "metallic" 
d-wave Bose liquid state instead of Bose condensing at a 
finite set of momenta. Such a phase is not accessible in a 
mean-field BCS analysis; however, it can be argued to be 
a reasonable outcome by observing that when deriving 
an effective boson Hamiltonian within perturbation the- 
ory at strong coupling, a ring term identical to the one 
we consider in the pure boson context is generated for 
increased hopping anisotropy between the two fermion 
species .'^Sl It is precisely this ring term that drives o ur bo- 
son system into a "metallic" DBL phase (see Sec. IIIB 
and Ref. ITS)) . Thus, such a "Cooper pair Bose metal" 
(CPBM) may exist at strong hopping anisotropy and in- 
termediate to strong attractive interactions. Evidence for 
such a phase was in fact recently found in a genuine at- 
tractive Hubbard model on the two-leg ladder in Ref. [221 
and our results here on three and four legs in the pure 
boson context may warrant further future studies of the 
CPBM, both theoretically and experimentally. 

A final line of motivation to study the DBL involves 
thinking of the hard-core boson ring model [see Eq. (|3|] 
as the easy-plane limit of an SU(2) invariant spin- 1/2 
model with four-spin cyclic ring exchange (see, for ex- 
ample, Ref. 1501) ■ III such a model, at zero magnetic field 
(half- filling in the boson language), exact diagonaliza- 
tion (ED) indicates, among other things, the presence 
of a spin-nematic uniaxial magnet intervening between 
a four-sublattice biaxial Neel state and a fully gapped 
SU(2) symmetric valence bond solid.l^ The same model 
has also been studied on the two-leg ladder.!^ Using the 
same multi-leg ladder approach we employ here, it would 
be interesting to study this model in the presence of a 
Zeeman field, which would correspond in the boson lan- 
guage to a finite chemical potential. Because the 2D 
DBL is expected to be generically stable only at densi- 
ties slightly below half-filling (and likely not present at 
exactly half-filling) it is conceivable that such a model 
in a sector of non-zero net magnetization could enter a 
spin liquid phase related to our DBL. Thus, understand- 
ing the physics of the DBL and the U(l) symmetric ring- 
exchange model of Eq. ([s]) represents a first step in un- 
derstanding this putative spin liquid phase. 

The rest of the paper is organized as follows. Section[ll] 
describes the basic machinery of our work, including the 
gauge theory description (Sec. HA) and the construc- 



tion of the variational wave functions (Sec. HB); this 
section also establishes our microscopic model (Sec. HC) 



and provides definitions of the physical measurements 
we consider (Sec. II D I. Section [TTT| addresses in detail 



the results of our study on the four-leg ladder, including 
most prominently the detection and characterization of 
the DBL [4, 2] gaplcss Bose metal phase. Section IV gives 
a similar analysis of our three-leg study, giving special 
attention to the DBL[3, 0] gapless Mott insulator phase. 
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This phase was the focus of our recent Letter, Ref.[Tni We 
first summarize the results of Ref. [19] and then present 
new, additional evidence and arguments in support of 
the stability of the gapless Mott insulating phase. Fi- 
nally, in Sec. [V] we discuss our conclusions and some of 
our plans for future work and extensions of the ideas pre- 
sented in this paper. Appendix \K\ presents our bosonized 
solutions to the gauge theories for both the DEL [4, 2] and 
the DBL[3, 0] as well as a stability analysis taking into ac- 
count the effects of short-range interactions. Appendix [B| 
contains results on the three-leg ladder at incommensu- 
rate densities and the successes and failures of the DEL 
framework in these systems. Specifically, we find strong 
evidence for a DEL[3, 1] phase for densities p < 1/3, while 
for p > 1/3 no DEL phase appears to exist; instead we 
have identified a phase consistent with a three-leg de- 
scendent of the "bond-chiral superfluid" phase predicted 
in a recent spin- wave analysis of the 2D J-K model.'SS In 
Appendix [C] we discuss the situation on the two-leg lad- 
der at half-filling and, in particular, why we fail to find 
a phase analogous to the DEL[3,0] in that system. 



II. PRELIMINARIES 

A. Gauge theory description: the DBL on the 
A^-Ieg ladder 

What follows is a generalization of what has been done 
in the appendix of Ref. 18 for the two-leg ladder. Here we 
will summarize the approach and state the key results. 
More details can be found in Appendix [A} 

We can describe the d-wave Bose liquid state by first re- 
expressing the bosonic operators as products of fermionic 
partons as in Eq. On the TV- leg ladder, each of the 
partons has the freedom to fill at most N ID bands in 
momentum space corresponding to the N transverse mo- 
menta; e.g., choosing periodic boundary conditions in the 
y direction, we have ky = 2j-K/N for j — 0, . . . , iV — 1. 
From our mean-field understanding of the state, we ex- 
pect the partons to occupy contiguous strips centered 
about fca; = in each band; such strips can best be 
described as "Fermi segments" and the edges of each 
segment as "Fermi points," the locations of which are 
the Fermi wave vectors. These points are the locations 
where partons can gaplessly be added or removed from 
the system and lead naturally to the concept of right- 
and left-movers (see, for example. Fig. [T]). Labeling the 

right-movers' momenta as k''p^\ where a = 1,2 refers 
to the two flavors of partons and ky labels the band, we 
can take the continuum limit in the longitudinal direction 
(i) and then linearize the dispersions around each Fermi 
point, capturing all the relevant low energy physics, and 
approximately decompose the partons as follows: 



dlix,y) ^ 
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iP[kpJx + ky{y-l) 



j(fey)t 



(6) 




37r/4 h 



-37r/4 [ 



FIG. 1: (color online), (top panel) Band filling for the 
DBL[4,2] showing the Fermi wave vectors for the right mov- 
ing di and d2 partons. (bottom panel) An overhead view of 
the occupied momentum states for the DBL[4, 2], highlighting 
the 2D nature of the phase; it is clear that the di Fermi sea 
is compressed along x, while the d2 Fermi sea is compressed 
along y. 



where P = R/L = +/— and denotes the right- and left- 
movers. Note that it is not necessary, in general, for 
either parton to occupy all TV bands; therefore, ky in the 
sum runs only over the partially filled bands for each 
flavor a. If the di partons partially occupy n bands and 
the d2 partons partially occupy m bands, we denote the 
resulting state as DEL[n,TO]. Since the two flavors of 
partons are interchangeable, we can always choose N > 
n > m. With this convention, the di (^2) partons are 
most easily associated with those that hop preferentially 
in the x (y) direction as described in the introduction. 
Finally, if we sum up the Fermi wave vectors for each 
parton, we recover the boson density p: 



iky] 

F2 



Nnp. 



(7) 



Continuing with our gauge theory description, we as- 
sign equal and opposite gauge charges to the di and ^2 
partons and turn on an appropriate U(l) gauge field; in 
quasi- ID, this behaves like a conventional electric field 
at long distances attracting the oppositely charged par- 
tons toward one another. In the limit of strong coupling, 
the partons are bound together on each site such that 
d\{r)di{r) — d\{r)d2{j^) = 6^(r)6(r), thus realizing the 
physical bosons. While the mean-field treatment would 
predict that the DEL[n, m] possesses n+m gapless modes 
(since the partons are completely free in this case), the 
effect of the gauge field is to render massive the overall 
charge mode (see Appendix |A]) reducing this number by 
one. Hence, a critical feature of the DBL[n, m] is the 
existence of n -f m — 1 gapless ID modes. 
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The gauge theory also offers an explanation for the 
singular behavior observed at incommensurate wave vec- 
tors in momentum-space correlators such as the boson 
momentum distribution function, rib(q), and the density- 



density structure factor, Db{q), for the DBL (see Sec. II D 



for explicit definitions of these quantities). Consider first 
nt,(q), the Fourier transform of the boson Green's func- 
tion, Gh(r). In the mean- field treatment, G^^(r) — 
Gf^ir)Gf/ir)/p; that is, simply the product of the in- 
dividual fermionic parton Green's functions. Therefore, 

will oscillate at various (Pk'p^^ + P'kp'^\Pky + 
P'k'y) wave vectors and decay as these momenta 

correspond to the gap less addition/removal of a boson, 
which involves creating/destroying a di and a d2 parton 
at certain Fermi points. When gauge field effects are 
taken into account, we expect these power laws to be al- 
tered for the different wave vectors with some oscillations 
dying off more quickly and some more slowly relative to 
the mean field. In Appendix we present a complete 
low-energy theory for the DBL phases on ladders, which, 
in principle, allows the exact calculation of all power law 
exponents. However, these exponents are also affected 
by short-range interactions and, due to the large number 
of parameters in the theory, we do not attempt to do 
these calculations explicitly. The gauge theory suggests 
a mechanism of Amperean enhancement whereby added 
(or removed) partons with opposite group velocities are 
favored since their opposite gauge charges generate paral- 
lel currents, which will subsequently attract one another 
thereby satisfying the gauge constraint and realizing the 
physical bosons. Assuming that the gauge interactions 
dominate in our theory, we can rely on the intuition of 
this Amperean rule despite the effects of the short-range 
interactions mentioned above. Thus we expect enhanced 
singularities in the momentum distribution function at 

wave vectors ±(k^pl^ — k'^p^ ,ky — k'y). As an example, 

the fermion bilinear rfiR'''^<i2L'^^ corresponds to such an 
enhancement. These "Bose points" are the quasi-lD fin- 
gerprints of the "Bose surfaces" of the parent 2D DBL 
phase; that is, the finite transverse momenta of the lad- 
der slice through the ID singular curves present in full 
2D resulting in these points of singular behavior. It can 
be shown within the context of the gauge theory that 



the other wave vectors, zt{k 



Fl 



^F2 



, h,y 



k'y), are al- 



ways suppressed relative to the mean field. In momen- 
tum space, nf,(q) shows clearly those wave vectors that 
are enhanced with sharp peaks or kinks in the curve. 

Turning now to the density-density correlator in real- 
space, this function can be approximated in the mean- 
field treatment as the average of the individual parton 
correlators: 



cr(r) 



Df^{v)]/2, and is 



therefore expected to oscillate at various "2/ci?" wave 
vectors and decay as l/x^. These can be thought of as 
particle-hole excitations between bands of the same fla- 
vor of parton. When the gauge field is turned on, the 
power law exponents will be modified and the Amperean 
rule predicts that excitations corresponding to aligned 



gauge currents will again be enhanced. Hence, we ex- 
pect enhanced singularities in the density-density struc- 
ture factor at wave vectors ±{k^p^ + kp^ , ky + k'y). As 

an example, the fermion bilinear c^^^^^^rfl^L"' corresponds 
to such an enhancement. The gauge theory predicts that 
the non-oscillatory, zero-momentum term in the boson 
density-density correlator will be unaffected by the gauge 
interactions, i.e., this term will simply remain and 

the other "2^^-" wave vectors, ±(fcp^^ — k^p^ ,ky — k'y), 
will always be suppressed relative to the mean field. As 
in the case of the boson Green's function, considering 
the momentum space density-density structure factor, 
-Df,(q), allows for a more organized view of those wave 
vectors that are enhanced by the presence of the gauge 
field. 

The gauge theory predictions involving the number of 
gapless modes and the locations of singularities allow us 
to identify the DBL[n, m] if and when it shows up as the 
ground state of a given model. We have proposed such 
a model in E q. ([s] ) and will discuss it in greater detail 
below in Sec. IHCT 

By solving the gauge theory model using bosonization, 
one can recast the DBL[n, to] as an interacting Luttinger 
liquid with rt + to — 1 ID modes with potentially distinct 
dispersion velocities and corresponding, nontrivial Lut- 
tinger parameters characterizing all power law exponents. 
This solution, along with a stability analysis, appears in 
Appendix [A] 



B. Wave function: DBL states on the A'^-leg ladder 

While the DMRG method can readily be employed on 
any quasi- ID model to extract out measurements of the 
ground state correlation functions, at best, one would be 
able to say that these functions decay as power laws and 
oscillate at incommensurate wave vectors if the ground 
state were indeed some DBL [n, to]. To make sense of 
the locations of these wave vectors, we have developed 
a more direct test for the presence of this phase using 
variational wave functions and VMC. These wave func- 
tions can be thought of as crude representations of the 
DBL[n,TO] phase as described by the bosonized parton 
gauge theory in the previous section. Following the mo- 
tivation of our mean-field parton description, we start by 
writing the bosonic wave function as a product of two 
Slater determinants wherein the orbitals occupied by the 
partons correspond to the band filling prescription laid 
out in the previous section (see Fig. [T]). But such a wave 
function without any further constraints would simply 
correspond to that of two fiavors of noninteracting spin- 
less fermions. We must project this wave function into 
the space of the physical bosons; i.e., the set of positions 
of the partons must be equivalent in both determinants 
[see Eq. ([T])] . This Gutzwiller projection of the wave func- 
tion plays the role of the strongly-coupled gauge field in 
the gauge theory. This sort of variational approach has 



FIG. 2: (color online). Schematic picture of how the anisotropy of the Fermi seas in the parton representation is driven by the 
ring-exchange term in the model, Eq. (111. If one allows the di and d2 partons to virtually hop along their preferred directions, 
an effective ring exchange for the bosons is generated. Because we are considering partons with fermionic statistics, the sign 
of the ri ng-exc hange interaction is rendered positive [cf. Eq. (13l]; in contrast, bosonic partons would generate negative ring 
exchange .'^^'^ The resulting sign can intuitively be understood by noting that the ring-exchange process involves one crossing 
of the partons, as depicted in the second leftmost plaquette shown above, so that the ring-exchange energy is minimized by 
having a positive (negative) coefficient for fermionic (bosonic) partons. 



been rather successful in ID systems, such as for the orig- 
inal studies of the t-J model.SsMI] 

The variational parameters in this wave function are 
simply how one defines the shape of the Fermi seas, built 
up from ID segments for di and d2 (see Fig. [T]). To be 
truly "d-wave" in nature, there must be some anisotropy 
in this filling such as in the DBL[2, 1], which was identi- 
fied and characterized in Ref. UHl or the DBL[4, 2] and 
DBL[3,0] states described later in this paper and in 
Ref. [TH But there exists much fiexibility in how one does 
this and many nontrivial fillings are possible in general, 
each corresponding to different locations of the singular 
features in the momentum-space correlators. Indeed, it 
is the agreement of these locations between the DMRG 
results and those of the most energetically competitive 
VMC states that forms the most compelling evidence for 
the existence of the DBL[n,m] in a given model. 

Additional variational freedom can be added by scaling 
the magnitude of the determinants in the boson wave 
function; this is one of the fruitful methods borrowed 
from the ID t-J studies mentioned above. Explicitly, if 
we start by writing schematically 



*6 ^ Vg [(detZ?i) X (detDa)] 



(8) 



where Di and D2 are the two matrices populated with 
e'^i ""j, the set {q|"^} is the Fermi sea for da, and Vg 
denotes the Gutzwiller projection, then we can replace 
each determinant as follows: 



detZ^a = \detDa 



detDa 
\detDa\ 



such that 



Vc 



detDi 



detD2 



IdetD] 



|i-pi 



IdctA 



|1-P2 



(10) 



The exponents pi and p2 can now be varied; this allows 
one to tune the m -\- n — 1 Luttinger parameters men- 
tioned in the previous section. We expect that the "bare" 
Gutzwiller wave function, i.e., Pi = P2 = Ij corresponds 
to all Luttinger parameters fixed at their trivial values 



of unity. Varying these exponents can also be seen as 
augmenting the pure Gutzwiller-projected Slater deter- 
minant wave function with a logarithmic Jastrow-type 
factor. 

Under the Gutzwiller projection, different mean-field 
states can lead to the same physical wave function (gauge 
redundancy). For our system and our practical purposes, 
this can be summed up as follows: if one shifts the Fermi 
sea of di by momentum Q and simultaneously shifts the 
Fermi sea of ^2 by momentum — Q, the wave function 
remains unchanged. The most significant consequence 
of this property is that there exists some arbitrariness 
in the choice of parton boundary conditions. Since we 
are going to be considering a finite system with periodic 
boundary conditions in both directions for the physical 
bosons, we actually have the option of choosing the par- 
ton boundary conditions, in each direction separately, to 
be both periodic or both antiperiodic. Either selection is 
mappable to the other by virtue of this shifting property, 
and so we will often make this choice based on what is 
most conceptually clear or aesthetically appealing. 



C. Microscopic model: hard-core bosons with 
frustrating ring exchange 

In the introduction, we stated the model of primary 
interest in Eqs. ([3])-([5]). The unfrustrated version of this 
model {K < 0), with both hard-core and soft-core bosons 
at and away from half-filling, has a long history in the 
past decade as a proposed candidate for harboring a de- 
confined quantum critical point and/or an exotic quan- 
tum phase dubbed the "exciton Bose liquid" (EBL),'^a 
relative of our DBL. However, for this particular model, 
both of these scenarios have been largely ruled out in a 
sequence of quantum Monte Carlo studies ,'^3'^ although 
recent work has shown that the EBL can be stabilized 
if one supplements a if -only model, Eq. dsj), with ring 
exchange on 1 x 2 and 2x1 plaquettes.'^"' 

Our focus here is on the explicitly frustrated {K > 0) 
case. This ring-exchange term and its role in manifesting 
the proposed DBL phase has a simple, intuitive physical 



7 



explanation, which is visuaHzed in Fig. [2] 

We can modify Eq. ([3| somewhat to allow for 
anisotropic hopping of the bosons; that is, 

H = Hhop + Hk, (11) 
Hbop = - J^(6|;6r+i + H.c.) - J_L J2iblbr+v + H.C.), 

r r 

(12) 

Hk =Kj2{blbr+ibl^,+ybr+y + U.c.). (13) 

r 

Now there are two dimensionless parameters, K/ J and 
J±/Jj in addition to the boson density p. Thus, for 
a fixed value of the density, we can explore a two- 
dimensional phase diagram in search of DBL phases. 

Our model lives on an A'^-leg square ladder wherein 
we can define lattice coordinates: x = 1, . . . , L^,, where 
Lx is the length of the chains, and y — 1, . . . , Ly, where 
Ly = N, the number of legs. The ring term in the Hamil- 
tonian applies to all elementary square plaquettes on the 
ladder. We use the convention in all numerical work that 
each hopping term for a given pair of sites and each ring 
term for a given plaquette is counted precisely once in 
the sums. 



D. Measurements 

We examine many quantities in this work and will de- 
fine them explicitly in this section. The first and most 
obvious is the ground state energy. This can be probed 
directly with DMRG and ED (for small system sizes). 
With VMC, we instead compute the trial energy of the 
variational wave function using our model [see Eq. (11)]. 



Minimizing this quantity with respect to the variational 
parameters (the band fillings and the exponents on the 
determinants), allows us to find the most competitive 
VMC state at each value of K/J and J±/J and plot 
out a VMC phase diagram, which can then be compared 
against DMRG results (see Sees. |III A]and|IV A[). 



We examine three main correlators to characterize the 
ground state for a given set of model parameters. First 
is the single-particle Green's function: 



(14) 



and its Fourier transform, the boson momentum distri- 
bution function: 

J-'xJ-'y , 

^ r If' 



where Ly is the number of legs of the ladder (interchange- 
able with iV in this paper) and is the number of sites 
in each chain of the ladder. 

Next, we define the boson density-density correlator: 



where p(r) = h'^(y:)h[v) and p = N),/ [L^Ly) (with no 
argument) is the average boson density. Also, its Fourier 
transform, the density-density structure factor, is defined 

as 



DM) 



1 



LxLy 

^ r.r' 



i?b(r,r')e' 



/^p'q•(I■-r') 



. (^^) 

We also found it useful to consider in some circum- 
stances current-current correlations on our ladder sys- 
tems. To this end, we define the current operator as 

Jf(r) = I {h\v + /i)6(r) - 6t(r)6(r + p)\ , (18) 

where p = x,y, and the current-current correlator as 



jf(r)Jnr': 



(19) 



The associated structure factor is simply the Fourier 
transform: 



Cr(q)^-^^Cr(r,r')e^'>- 



(r-r') 



(20) 



In this work, we have only considered the y-y and x-x 
correlations, i.e., Cj'^(q) and C^'^(q). 

As mentioned earlier, having a means by which to de- 
termine the number of gapless ID modes is a critical di- 
agnostic for detecting a DBL[n, m] state. We can refer to 
conformal field theory (CFT), the results of which sug- 
gest a relationship between the number of gapless modes 
and the central charge, a quantity that can be extracted 
from the scaling form of the entanglement entropy. For 
a subsystem size AT in a system of overall length with 
periodic boundary conditions, the von Neumann entan- 
glement entropy S is given by the scaling fornP^ 



S{X,Lx) = -log 



Lx ttA 
— sm — — 



(21) 



A(r,r')^(b(r)-p][p(r')-p]) 



(16) 



where A is a nonuniversal constant independent of the 
subsystem length and c is the central charge that we seek. 
In our ladder studies, we only consider clean vertical cuts 
through all N = Ly chains so that the left block contains 
NX sites {X out of L^ total rungs). Despite the fact that 
the bosonization analysis leads to a theory of ti -f to — 1 
free bosonic modes with generally different velocities such 
that the full system is not conformally invariant, we still 
expect the overall measure of gaplessness to be insen- 
sitive to this detail: since we are extracting c from the 
ground state wave function only, which has no knowledge 
of mode velocities, the scaling form of S should not de- 
pend on the velocities. This is consistent with the known 
"L log L" scaling for free fermions"* ^ and projected Fermi 
sea stateJi^ in 2D, in which conformal invariance is not 
present. By measuring S using DMRG while varying A, 
we can attempt to fit this form to the data and extract 
the constant fit parameters c and A. Doing so requires 
highly converged DMRG data on large system sizes, L^, 
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in order to get reliable estimates of the central charge. 
Since this task becomes increasingly more computation- 
ally challenging with greater spatial entanglement, our 
DBL [n,m] phases can be particularly difficult to analyze 
as c — n-\-m—l grows large. For example, while we were 
successful in measuring a central charge of approximately 
two for the DBL[3, 0] (see Sec. IVB), we exhausted the 



computational resources for the DBL [4, 2] (see Sec. Ill B ) 
where we expect c = 5. 

Finally, we also considered comparisons of the ground 
state momentum at various points in the phase diagram 
for the DBL [3, 0]. We looked at small system sizes where 
finite-size effects are expected to be significant and ex- 
tracted the total momentum from ED calculations, which 
is easily inferred from ED results since the diagonaliza- 
tion is carried out in blocks diagonal in the momentum 
quantum number. We then compared these results to 
the ground state momenta of the VMC states, which are 
also trivial to compute since we are using complex- valued 
orbitals, e''^ '': one simply sums up all of the individual 
momenta of the partons, which fill definite momentum 
orbitals (see Figs. [5]and|6j bottom panel). 

The DMRG calculations presented in this paper were 
performed with anywhere between D ~ 1000-9000 states 
per block and at least 6 finite-size sweeps, where each 
"sweep" traverses the L^-Ly sites of the lattice twice. 
The accuracy of the results is strongly dependent on the 
phase being studied, the chosen system size and bound- 
ary conditions, as well as the physical quantities being 
measured. For example, within the four-leg DBL[4, 2] 
phase discussed in Sec. |IIIB| the momentum space corre- 
lators on which we focus, nh(q) and I?b(q), are converged 
to a relative error of less than 10~^, and the ground state 
energy to a relative error on the order of 10~*; the den- 
sity matrix truncation error is on the order of 10"^. In all 
other identified phases, the accuracy of our results is bet- 
ter; e.g., within the four-leg superfluid phase (see Fig. [4]), 
the truncation error is on the order of 10~* keeping only 
D = 2 000 s tates, and in the thrcc-lcg DBL[3,0] phase 
of Sec. llVBl the truncation error is on the order of 10^^ 
keeping D = 4000 states. These quoted errors are for the 
case of fully periodic boundary conditions, which we pre- 
fer due to the incommensurate nature of the DBL phases. 
The entanglement entropy is typically the last quantity 
to converge, and obtaining highly converged entropy data 
in the multi-mode critical systems we encounter is an ex- 
tremely challenging numerical task. We were able to ob- 
tain such converged entropy data (used for determination 
of the central charge, c) within the superfluid, DBL[3, 0], 
and DBL[3,1] phases, but were unfortunately unable to 
do so within the DBL[4, 2] phase. Finally, we note that 
in such (potentially) "very critical" quasi-lD phases, it is 
not possible to go to very long systems (L^, > 100, say, at 
fixed Ly = N ^ 3,4) with the DMRG to definitively rule 
out eventual small gaps and corresponding long (finite) 
correlations lengths; this is further compounded by the 
fact that the finite bond dimension matrix product states 
produced by DMRG give exponentially decaying correla- 



tions by construction. However, such weak instabilities 
should not be expected a priori, and for the DBL[4, 2] 
and DBL[3,0] phases, we have done the best we can to 
rule out such scenarios by going to reasonably long sys- 
tems using both fully periodic and cylindrical boundary 
conditions. 

At this point, it is convenient to mention a subtlety 
regarding the DMRG and how it handles ground states 
with non-zero momentum. To begin with, it is not possi- 
ble to directly extract the ground state momentum from 
the DMRG states since the real-space blocking construc- 
tion necessarily breaks translational invariance. If this 
momentum is indeed zero or an integer multiple of tt, 
there is no ambiguity in our DMRG wave function and 
the choice of real- versus complex-valued wave functions 
is irrelevant. If, however, this is not the case and a given 
ground state has momentum Q, then it necessarily has a 
time-reversed partner with momentum — Q. The DMRG 
ground state is thus some real-valued combination of 
these two states. While the lattice-space measurements 
discussed above may depend slightly on the details of this 
combination, the differences become less significant with 
increased system size. In all the VMC/DMRG compar- 
isons below, we choose points in the phase diagram where 
Q = so as to completely avoid this ambiguity. 

With the preliminary details having now been fully 
discussed, we are now prepared to give detailed descrip- 
tions of the results for the model, Eq. (Ill, on three- and 
four-leg ladders. 



III. GAPLESS BOSE METAL PHASE ON THE 
FOUR-LEG LADDER 

Here we present results of our VMC and DMRG study 
of the model, Eq. (11), on the four-leg ladder {N = 
Ly — 4) with periodic boundary conditions in both the 
transverse and longitudinal directions. Naturally, we first 
searched for metallic DBL[n,m] phases on the three-leg 
ladder but were unsuccessful in finding any that were ex- 
tensible to 2D, most likely due to the non-bipartiteness 
of the lattice. We did find some other interesting results. 



which are discussed in Sec. |IV] and Appendix [B] Here, 
we focus on X 4 system sizes with = 12, 18, 24. 
We expect qualitatively different DBL[?t,, m] states in the 
two naturally defined density regimes: < p < 1/4 and 
1/4 < /9 < 1/2. For small boson densities in the former 
case, we find that the system readily phase separates as 
the ring coupling is increased. For this reason, and that 
we would like to find metallic DBL[n, m] phases that are 
extensible to 2D, we considered densities in the latter 
regime. While the behavior is expected to be generic for 
p between 1/4 and 1/2, for concreteness we will use boson 
density p — 5/12 for all three system sizes. The results 
are summarized in the phase diagram. Fig. (3) which we 
shall discuss first. The differently colored regions are de- 
termined from a VMC study on the 18 x 4 system while 
the DMRG points are from all three system sizes. Fol- 



9 



3 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

OSF 

- ° [4,2] oooflnnoDfloo- 
♦ PS 

ooonnnnnfl«« 

SF 

□ [4,2] 

-OOBnnDDDI^««0 



01234567 
K/J 

FIG. 3: (color online). Phase diagram of the four-leg system 
at boson density p — 5/12, using system sizes 12 x 4 with 
A^i, — 20 bosons, 18 x 4 with Nt = 30 bosons, and (for some 
points) 24 X 4 with Nb — 40 bosons. The colored regions are 
delineated using VMC data. The gray region indicates phase 
separation and is delineated schematically by considering the 
boson density in real space with the DMRG approach. DMRG 
points for the superfluid, DBL[4, 2], and phase separation are 
indicated by white circles, green squares, and gray diamonds, 
respectively. Finally, in the region labeled "[4,4]," the VMC 
finds a DBL state with four equally occupied bands for both 
the di and d2 partons, although this region is likely just a 
superfluid. 



lowing that, we shall give an in depth analysis of the 
DBL [4, 2] phase. 



A. The four-leg phase diagram at p — 5/12 

For small K/J, the DMRG confirms the existence of 
a quasi- ID version of the generic superfluid described 
earlier (points marked with circles in Fig. |3]) exhibiting 
quasi-long range order with the bosons "condensed" at 
q = (see Fig. 4j top panel) and a central charge of 
c ~ 1 corresponding to the single gapless mode, as de- 
termined by measuring the entanglement entropy. We 
model this phase in the VMC calculations using a simple 
Jastrow wave function that simulates the inter-particle 
repulsion using a constant potential for nearest neigh- 
bors and a power law in the separation distance in all 
other cases. The Jastrow wave function contains three 
floating point variational parameters that are optimized 
by minimizing the trial energy over a finely grained mesh 
of the phase diagram. The phase diagram indicates that 
this modeling is largely successful as it reproduces rea- 
sonably well the phase boundary to the DBL [4, 2] regime. 
The density-density structure factor for a typical point 
in the superfluid phase is shown in Fig. [4j bottom panel. 
The characteristic behavior around qx = for qy = 
is clearly present. 

Along the cut J± = J, the DBL [4, 2] onsets for roughly 
K/J > 1.75 via a first-order phase transition. This is de- 
termined using DMRG (points marked with squares in 
Fig. |3]) by looking for the appearance of singular features 
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FIG. 4: (color online). Boson momentum distribution func- 
tion (top panel) and density-density structure factor (bottom 
panel) for the superfluid point at J±/J — 1, K/J = 1 for 
a 24 X 4 system size with Nb — 40 {p — 5/12), as obtained 
by DMRG. The values of ni,(q) at Qy — ±n/2,n have been 
scaled up by a factor of 20. The q = condensation is readily 
apparent in n(,(q), as is the \qx\ behavior of -Di,(q) near = 

at Qy = 0. 



at incommensurate wave vectors in the boson momen- 
tum distribution function and the density- density struc- 
ture factor (see Figs, [s] and [6] for characteristic points). 
The precise locations of these features evolve as K/ J and 
J±/J are varied within the DBL region of the phase dia- 
gram (see Figs. [5] and [6|. This is consistent with the par- 
ton picture where larger K/ J leads to greater anisotropy 
in the Fermi seas and hence the number of c?i partons in 
the ky band N^'^^'^ Ni^/A for all ky, which drives the 
evolution of the peak locations. In the extremal state, 



^dl^' — which can only be truly realized when Nj, 

is divisible by four, n{,{q) becomes independent of qy in 
the variational wave function since there is a conserved 
number of bosons in each chain (see Fig.l6|. With VMC, 
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we use the wave functions as described in Sec. IIIBI to 
model the DBL. At the onset, we perform an exhaustive 
search over all possible, unique band filling configura- 
tions, the only constraint being that the filled orbitals 
form a contiguous strip centered at fc^; = within each 
band. Note that there are cases where the band occu- 
pation numbers are such that no choice of parton lon- 
gitudinal boundary conditions results in all bands being 
symmetrically filled. In such cases, the resulting wave 
function has a non-zero total momentum in the x direc- 
tion and all possible, unique resolutions of the "leftover" 
particles are considered in the search. This first analy- 
sis is performed with the "bare" wave function, i.e., with 
the exponents on the determinants, pi and p2, set equal 
to one. Using these results, a pool of the most competi- 
tive bare VMC states within the parameter regime shown 
in Fig. [3] is formed and then a second analysis stage is 
performed varying these exponents over < {^1,^2} < 1- 
The best DBL state is compared energetically to the best 
Jastrow state at each point in a finely grained mesh over 
the phase diagram and this is how the VMC phase bound- 
ary is determined. 

For still larger values of the ring coupling, roughly 
K/J > 4.25, the DMRG data suggest spatial phase 
separation into regions of zero density and density 1/2 
(points marked with diamonds in Fig. |3|. Specifically, 
our DMRG generally "gets stuck" in a nonuniform state 
with the (1 — 2p)Lx centermost rungs at near zero den- 
sity and the remaining 2pLx outermost rungs at 1/2 den- 
sity; there is a corresponding sharp feature in the density- 
density structure factor, Db(q), at q = ±{27^ / Lx^O). The 
actual ground state is presumably an equal superposition 
of such states with the center of the low-density region 
at all different rungs of the ladder, to which we could 
never hope to equilibrate with the DMRG. This tendency 
to increase the local density is expected since the ring 
term induces an effective attraction of the bosons. It 
is apparently the case that the hopping terms are suffi- 
cient to stabilize phases, such as the superfluid and the 
DBL [4, 2], over a large region of the phase diagram for 
smaller values of K. 

Throughout our phase diagrams (see Fig. [s] above and 
Figs. [7] and 14 below), question marks denote a lack of 
surety in the DMRG data for specifying a phase. In 
most cases, however, these points are located near phase 
boundaries on which we place limited focus in this paper. 



B. The DBL[4,2] phase 

We now turn our attention to the DBL[4, 2] phase it- 
self (see Fig. [ij, which is a conducting, metal-like phase 
with one more gapless ID mode than the number of legs 
on the ladder. This phase shows no signs of ordering and 
breaks no symmetries; in particular, it respects the inver- 
sion symmetry of the lattice. We choose two character- 
istic points deep within this region of the phase diagram 
for further analysis: J±/J = 1,0.1 at fixed K/J = 2.5. 



We scale up the most competitive VMC states at these 
points from the 12 x 4 system by multiplying all band 
occupation numbers by two in the hopes of accurately 
capturing the behavior of the phase on a 24 x 4 system 
with Nb — 40. The same points on this larger system size 
are studied with DMRG and the comparison of these re- 
sults is shown in Figs. [5]and[6j We perform this scaling 
procedure to avoid doing a full VMC energetics study 
on the 24 X 4 system where the pool of potential trial 
wave functions is sufficiently large as to render the study 
intractable. Naturally, the VMC method is capable of 
calculating correlators for a single state on system sizes 
much greater than this, while the DMRG approach is 
reaching its computational limit. 

We now describe in detail the state found at J±/ J = 1 
and K/ J = 2.5 (see Fig. [5| as a representative of the 
DBL [4, 2] phase and give an explanation for the loca- 
tions of the singularities in the structure factors consis- 
tent with the gauge theory predictions. First, we will 
take antiperiodic boundary conditions for both partons 
in both the x and ij directions so that the bands will 
be filled symmetrically; we remind that this is consistent 
with the physical bosons obeying periodic boundary con- 
ditions in both directions. The band filling situation is 
displayed visually in Fig. [5j bottom panel; clearly, this 
state has zero total momentum. For the c?i partons, the 
ky — ±7r/4 bands each have 12 orbitals filled, while the 
ky = ±37r/4 bands each have 8 orbitals filled. This leads 



to Fermi wave vectors of k 



(±7r/4) 
Fl 



127r/24 = 7r/2 and 



^(±37r/4) ^ 8,^/24 = 7r/3. (We point out that there is 
a technicality here: in order to be consistent with the 
continuum limit, it is best to define the Fermi wave 
vectors as halfway between the last filled orbital and 
the first unoccupied orbital.) For the d2 partons, only 
the ky — ±tt/4: bands are occupied and have 20 par- 
ticles each; hence the relevant Fermi wave vectors are 

k'^P2^^^ — 207r/24 = 57r/6. One can verify explicitly that 
the sum rule of Eq. ([?]) is satisfied for both parton fla- 
vors. The optimal variational exponents on the determi- 
nants for the corresponding point in the 12 x 4 system 
are pi — P2 = 0.8, which are close to the bare values 
of unity; these exponents have been carried over to the 
scaled up 24 x 4 state. 

Looking first to the boson momentum distribution 
function (Fig. [sj top panel), we see that two of the 

predicted enhanced momenta show up: ±(A;^'^^^'' — 



4±//^\0) = ±(^/3,0), ±(4t^'' 
±(77/2, tt). The other two that are predicted to be en- 
hanced are present in the VMC data, although weak, 
and appear to be smoothed out in the DMRG data: 
±(4t^^)-4=FV4)^^^/2) = ±(7r/3,±7r/2), ±(41^/'^" 
^(±37r/4)^-p^^2) = ±(7r/2,=F7r/2). Note that due to the 
inversion symmetry of the lattice, the ky — ±7r/2 bands 
are always degenerate in all cases. While the amplitudes 
of the peaks in the VMC data are not quite right, the 
agreement of the singular locations is striking. The fact 
that the coincident singularities at qy — Q^tt are quite 
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FIG. 5: (color online). Boson momentum distribution func- 
tion (top panel) and density-density structure factor (mid- 
dle panel) for a characteristic DBL[4, 2] point: J_l/J = 1, 
K/J = 2.5. The system size is 24 X 4 with iVft = 40 (p = 5/12). 
DMRG data are joined with solid curves while the VMG data 
are joined with dashed curves, (bottom panel) Schematic rep- 
resentation of the energy-optimized VMC state for this point. 
The circles represent momentum orbitals in each of the four 
bands; filled circles denote occupied orbitals. Each band is 
centered about = Q and we are using antiperiodic bound- 
ary conditions in both the x and y directions for the partons. 



FIG. 6: (color online). The same quantities are plotted as 
in Fig. |5] for a 24 X 4 system with A'^j, = 40 bosons (p = 
5/12), except now the measurements are calculated at the 
point Jx/J = 0.1, K/J = 2.5. The bottom panel depicts 
the new energy-optimized VMC state, whose di configuration 
differs from that of Fig. [5] Note the feature in the VMC data 
for ni,(q) at — ±37r/4, which is a "suppressed" singularity 
resulting from creating a boson on the same side of the di 
and d2 coincident bands: ±{kFi + kpT) = ^Sn/i. Analogous 
features are also present in the VMC data of Fig. [5] As well 
as demonstrating the ability to vary the location of the Bose 
points within the phase diagram, these results in conjunction 
with those of Fig. [5] also highlight the general nonuniversality 
of the amplitudes of the power law singularities. 
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strong in the DMRG, even more so than in the VMC, is 
itself encouraging for stabihty of the DBL phase. 

Considering next the density-density structure factor 
(Fig. [5] middle panel) , the qualitative and quantitative 
agreement between the VMC and DMRG data is excel- 
lent. In general, we do not expect the peak amplitudes of 
the VMC results to match those present in the DMRG 
data; these depend on the details of the gauge theory 
and the many input parameters about which the trial 
wave functions have no knowledge. This is the situation 
with the boson momentum distribution function data, 
whereas the -D6(q) curves are actually nearly coincident. 
However, the agreement of the locations of the peaks is 
what matters most in giving credence to the parton de- 
scription of the phase. The DMRG results pick up half 
of the predicted enhanced wave vectors, although some 
of these features are smoothed out in the VMC results: 
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tually shows up clearly in the DMRG but not in the 
VMC, where it is smoothed by inclusion of the varia- 
tional exponents [see Eq. (10)]. The VMC alone sees one 
of the other four wave vectors at which the gauge theory 
predicts enhanced singularities, but this does not show 

up in the DMRG data: ±(41''^''^ + T7r/2) = 

±(57r/6, =F7r/2). The remaining three wave vectors are 



not seen in either data set: ±(fc 
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In Fig. [6l we show boson momentum distribution 
and density-density structure factor measurements as ob- 
tained by DMRG and VMC at one more point in the 
phase diagram: J±/J = 0.1 and K/J = 2.5. Again, 
we take a 24 x 4 system with A'^;, = 40 bosons. Now, 
the di configuration consists of four equally filled bands, 
i.e., each band ky = ±7r/4, ±37r/4 contains iVb/4 = 10 
partons, while the d2 configuration is the same as that 
in Fig. [5j That is, our Fermi wave vectors are now 
given by kpi = 41"^*^ = 41^"^*^ = 10^/24 = 5^/12 
and kF2 = ^^F2'^^ = 207r/24 = 57r/6. This varia- 
tional wave function has a conserved number of bosons 
in each chain and hence nb(q) is independent of qy and 
Db{qx = 0,Qy) = for all g^; these properties hold only 
approximately in the DMRG due to small fluctuations of 
the single-chain boson number for J± 7^ 0. As before, the 
locations of the main singular features are consistent with 
what we expect from the gauge theory. Unlike in Fig. [5] 
we now see clear features in the DMRG data for nf,(q) at 
Qy = ±7r/2, which is to be expected since 7if,(q) is approx- 
imately independent of qy] however, the features in D{,(q) 
at qy = ±7r/2 are still smoothed in the DMRG data while 
being present, but very weak, in the VMC data. Note 
again the feature in -Db(q) at ±{2kpi,0) — ±(57r/6,0), 
which is predicted to be enhanced coming from the gauge 



theory and which actually shows up in the DMRG with- 
out being present in the VMC. Overall, the combined 
results of Figs. [5] and [6] clearly illustrate the tunability 
of the locations of the Bose points within the DBL [4, 2] 
phase. In going from Fig. [5] to Fig. [6] we have chosen to 
decrease J± instead of increasing K to induce evolution 
of the Bose points because (upon increasing K at fixed 
J±/J = 1) the system phase separates before entering 
the di configuration realized in Fig. [6] 

Considering Figs. [S] and [6] as a whole, the fact that 
not all predicted features are visible is not surprising nor 
does it detract from our argument in support of the re- 
alized phase being a DBL [4, 2]. In reality, one must con- 
sider all allowed interactions in the full gauge theory, and 
due to the complexity of this multi-mode Luttinger liq- 
uid, it is impossible to know how the short-ranged inter- 
actions will affect the anomalous power law exponents. 
The gauge theory, insofar as we can interpret it, only 
predicts potentially enhanced wave vectors, as well as 
those that are definitely suppressed, i.e., those that do 
not satisfy the Amperean rule. Also, it is important to 
note that the amplitudes of the power law singularities 
are nonuniversal quantities which we could never even 
hope to predict within our gauge theory framework. It 
is entirely possible, perhaps even likely given the over- 
all agreement between the VMC and DMRG, that such 
matrix element effects could be responsible for the ab- 
sence of some features in the DMRG momentum space 
correlators depicted in Figs. [5] and |6) note that even the 
amplitudes of the corresponding VMC singularities are 
relatively small in these cases. It has been suggested 
that similar effects may be playing a role in smoothen- 
ing some singular features in the four-leg spin Bose metal 
phase (so-called SBM-3) of Ref. HSl— note that DBL [4, 2] 
and SBM-3 are closely analogous phases realized in the 
respective ring models considered here and in Ref. I15[ 

We conclude this section by discussing scaling of the 
entanglement entropy within the DBL [4, 2] phase and our 
efforts to try to extract an effective central charge c = 5. 
Our gauge theory analysis of DBL [4, 2], as presented in 
Appendix |A 1[ predicts c = 4-|-2— 1 = 51D gapless 
modes, and in principle, we should be able to extract 
this by examining the scaling of the entanglement en- 
tropy according to Eq. (21). However, in practice, this is 



an exceedingly difficult task due to the large amount of 
spatial entanglement in the ground state as implied by 
such a phase. Specifically, to fully converge the von Neu- 
mann entanglement entropy the required D, which is the 
number of reduced density matrix basis states retained 
in the DMRG, becomes prohibitively large given current 
computational resources. 

We have only been able to attain full convergence on 
small system sizes; however, the results on these small 
systems are indeed suggestive of c ~ 5. For example, 
fitting entropy data of an 8 x 4 system with Ni, — 12 
bosons at the DBL[4,2] point Jj_/J = 1, K/J = 2.5 
gives c = 4.85 ± 0.05, where the error quoted is the error 
in the fit only (the entropy data itself is converged to 
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about 10~^ keeping D = 8000 states). Unfortunately, we 
have been unable to fully converge the entropy data for 
> 8 systems and have thus not been able to get reli- 
able estimates of c for reasonably long systems. For 12x4 
systems, the entropy [say at the middle of the system, 
S{X = La:/2, Lx)] is still generally slowly increasing even 
when keeping up to D = 9300 states, although fits to the 
still unconverged data give c > 4. The very slow increase 
in S with D is to be expected since S is known to grow 
only logarithmically with D for ID critical systems.'^ At 
some points in the phase diagram on the 12 x 4 system, 
the convergence of S is better due to a smaller value of 
the constant A in Eq. (21). In these cases, c is even 



closer to 5; e.g., at the point J±/J = 1.5, K/J — 2.5 
with Ni, = 20 bosons on a 12 x 4 system, we get c > 4.6 
with the entropy still very slowly increasing. We caution 
that in our experience with other realized DBL phases, 
e.g., DBL[2, 1] on two legs (see Ref. W), DBL[3, 1] on 
three legs (see Appendix B 1 1, and DBL [3, 0] on three legs 
(see Ref. [inland Sec. IV), small systems tend to overes- 
timate c from its expected value, usually on the order of 
about 10%. However, DBL[4, 2] is a qualitatively distinct 
phase from those mentioned above, so such trends may 
no longer apply here. Finally, we have tried using open 
boundary conditions (OBC) in the x direction (cylindri- 
cal boundary conditions overall) in which case it is much 
less costly to obtain converged entropy data on larger 
systems; e.g., we were able to get the entropy converged 
to 1% on systems up to 48 x 4 by keeping D — 6000 
states. However, this data is plagued by very strong os- 
cillations originating from the finite wave vectors clearly 
present in the DBL [4, 2] phase, and these oscillations ul- 
timately make an accurate determination of c impossible 
in this case. All in all, even though we have been un- 
able to definitively confirm c = 5 gapless modes with the 
DMRG measurements, our calculations in no way rule 
out the possibility of this result. 



C. Discussion 

It is clear that the agreement of the locations of sin- 
gularities in the structure factors is excellent, although 
there remain quantitative differences. It is important to 
bear in mind that the Gutzwiller wave functions provide 
only a caricature of the DBL[4, 2] phase (e.g., they should 
not be expected to accurately capture long-distance prop- 
erties of the underlying gauge theorji^H) and it is not the 
goal of this study to find the exact ground state of our 
ring model using a variational wave function; after all, we 
already have the ground state using DMRG. We again 
emphasize that the lack of some of the expected features 
in the DMRG measurements is not at all surprising. The 
full DBL [4, 2] theory, with both gauge and short-range in- 
teractions, is an extremely complicated multi-mode Lut- 
tinger liquid with five potentially nontrivial Luttingcr pa- 
rameters. Thus, the manner in which the power law ex- 
ponents is altered by the allowed short-range interactions 



is difficult to assess, as is the ultimate effect of nonuni- 
versal power law amplitudes. Overall, we find the data 
presented here to be compelling evidence for the validity 
of our parton description of the remarkable phase found 
in the four-leg frustrated J-K model and its correct iden- 
tification as a quasi- ID descendent of the d-wave Bose 
liquid. Furthermore, in Appendix |A 1[ we consider the 
effects of short-range interactions on the stability of the 
DBL [4, 2] within the scope of the bosonized gauge the- 
ory and conclude, from an analytical perspective, that 
DBL [4, 2] is likely a stable quantum phase. 

Further numerical confirmation could have in princi- 
ple been obtained by examining the entanglement en- 
tropy in the DMRG and extracting the expected cen- 
tral charge, c = 5 = 4-1-2 — 1. However, this study 
turned out to be inconclusive, mainly due to an inability 
to converge the entropy on large periodic systems, but 
also because open (cylindrical) systems are plagued by 
oscillations originating from the finite Bose wave vectors 
present in our DBL theory. Going forward, we believe 
that there is much room for improvement, both analyt- 
ically and numerically, using entanglement scaling prop- 
erties to characterize these systems. On the numerical 
front, it would be particularly interesting to use recently 
developed Monte Garlo techniques^® to compute the 
Renyi entropy (as opposed to the von Neumann entropy) 
directly in our projected variational wave functions; aside 
from being interesting in its own right, this would also 
give us a guide as to how much entanglement to expect in 
the DMRG calculations. For the Renyi entropy, we may 
expect subleading oscillatory corrections to scalin^^H for 
the DBL phases, again due to the presence of finite Bose 
wave vectors, and this could make the data more difficult 
to analyze. Finally, it may prove easier to access long- 
distance properties, e.g., the central charge, of nontrivial 
quasi-lD gapless phases such as our DBL[4, 2] and the 
proposed SBM-3 phase of Ref. [THI by using recentl y de- 
veloped "entanglement renormalization" techniques j'^^'^ 
an endeavor we leave for future work. 

We now compare and contrast the nature of the 
DBL[4, 2] to that of the DBL[2, 1] on the two-leg ladder 
characterized in Ref. 18, The former can be seen as sim- 
ply a scaled up version of the latter; the DBL[4,2] has 
the same band picture as the DBL[2, 1] but with pairs 
of coincident bands replacing the single bands (compare 
Fig. [l]of this paper to Fig. 2 of Ref. [15]). In this sense, 
the DBL[4, 2] is precisely the state we would expect in 
moving toward 2D. The most significant achievement of 
the work presented in this paper over what has been done 
previously is the realization of a DBL [n, to] phase where 
n+m—1 > N, that is, where the number of gapless modes 
exceeds the number of legs on the ladder. This situation 
has an interesting consequence: since n < A'^, to > 1, 
necessarily, and therefore the determinant for the d2 par- 
tons cannot be interpreted as a simple Jordan- Wigner 
string multiplying the di determinant. For m — 1, the 
c?2 determinant enforces a condition of no more than one 
particle per rung; in this case, the physics of this deter- 
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minant is essentially ID, allowing for the Jordan- Wigner 
interpretation. But when m > 1, the d2 determinant has 
a more subtle effect on the sign structure of the boson 
wave function. In Ref. [TSJ the authors searched for a 
stable DBL[2, 2] on the two-leg ladder, which also would 
have satisfied this condition. The bosonized gauge the- 
ory suggests a strong instability toward an s-wave paired 
phase, and the DMRG confirmed that this was indeed the 
case. Thus, no DBL[2, 2] was ever stabilized. Finally, we 
point out that the prevalence of phase separation in the 
phase diagram (see Fig. |4]) does not seem to have grown 
in going from A'^ = 2 to 4 legs, which is encouraging for 
the eventual stability of the DBL in 2D. 

Prior to studying = 4, we thoroughly explored the 
three-leg ladder in the hopes of finding a stable DBL [3, 2] 
phase, which would have also satisfied n~\-m—l > N , for 
densities greater than the commensurate value p > 1/3. 
With our chosen model, Eq. ( 11 1, we were not successful 



and instead found what appears to be a three-leg de- 
scendant of an unusual superfluid predicted in a recent 
spin-wave treatment of the J-K model on the 2D square 
lattice. •^■^ Although this phase can be qualitatively un- 
derstood with an exceedingly simple classical analysis 
of Eq. (Ill, it is itself rather exotic as it is a gapless 



phase that breaks both time reversal and translational 
symmetry (see Appendix |B 2| for details). We did, how- 
ever, find a stable DBL [3, IJ at density p = 1/4 < 1/3, 
which is a natural extension of the DBL [2, 1] on two legs 
and which we describe in greater detail in Appendix |B 1[ 
The DBL[3, 1] has exactly as many gapless modes as the 
number of legs, same as the DBL[2, 1]. The study of the 
three-leg ladder also led to our discovery of a novel, truly 
quasi- ID, insulating phase with two gapless modes and 
a d-wave sign structure in the wave function; that is, a 
DBL [3,0], which exists only for the special, commensu- 
rate case p = 1/3. We now discuss this phase further in 
the following section. 



paper self-contained. 



A. The three-leg phase diagram at p = 1/3 

The three- leg phase diagram at 1/3 filling is presented 
in Fig. [7] where we show VMC and DMRG results for a 
24 X 3 system size with Ni, = 24. There are three major 
phases plus some other regions where our understanding 
is currently limited. First, a large portion of the diagram, 
for small to moderate values of K/J and J±/J, is the 
generic superfluid phase, which has the same properties 
as in the four-leg case above. The DMRG points that 
we have identified as SF are marked with circles. Again, 
the DMRG confirms a central charge of c ~ 1 consistent 
with the predicted single gapless mode. This phase is 
modeled in the VMC using the same type of Jastrow 
wave functions as in the four-leg case. 

For larger values of K/J with J±/J < 2, we see evi- 
dence for the gapless Bose insulating phase: the DBL[3, 0] 
(marked with squares in Fig.[7j see Figs, [s] and |9] for char- 
acteristic points) . This phase is a particularly interesting 
incarnation of the DBL in that one parton (di) is gapless 
while the other parton (1^2) is gapped, the latter being 
due to a fully filled ky = band at p — 1/3. Due to the 
presence of this fully filled ^2 band and the absence of any 
partially filled bands (a situation that can only occur at 
commensurate densities; see Sec. IV B[ in particular the 
bottom panel of Fig. |9] for more details), there exists a 
gap to adding a boson and hence we expect the real-space 
single-boson Green's function to decay exponentially, a 
prediction that is confirmed by the featureless boson mo- 
mentum distribution function (top panel of Fig. [9]). The 
density-density structure factor (middle panel of Fig.|9|, 
however, shows singular peaks at incommensurate wave 
vectors; these are, in fact, the signatures for which we 



IV. GAPLESS BOSE INSULATOR PHASE ON 
THE THREE-LEG LADDER 

In this section, we present DMRG, VMC, and ED re- 
sults for our study of the same model [see Eq. (11)] on 
the three-leg ladder {N ~ Ly = 3) with periodic bound- 
ary conditions in both the longitudinal and transverse 
directions. Primarily, we focus on the DBL[3,0] phase, 
which is stable at boson density p — 1/3. This phase was 
characterized and discussed in detail in our recent Letter, 
Ref.TQ, in which it was referred to as a gapless Mott insu- 
lator (GMI). For completeness, we now summarize those 
results as well as present additional and compelling data 
and arguments in support of this exotic state as a stable, 
quasi- ID phase realized by our ring model. First, we will 
tour the phase diagram for a model finite-size system, 
and then we will focus the discussion on the DBL[3, 0]. 
Please note that Figs.[7]and|9]were presented in our previ- 
ous work,^^ but are reproduced here to make the present 
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FIG. 7: (color online). Phase diagram of the three- leg system 
at boson density p = 1/3, using a system of size 24 x 3 with 
A^i, = 24 bosons. The colored regions are delineated using 
VMC data; the white regions indicate where our understand- 
ing of the phase diagram is limited. DMRG points for the 
superfluid, rung Mott insulator, and DBL[3,0] are indicated 
by blue circles, green squares, and yellow triangles, respec- 
tively. This figure has been reproduced from Ref. (HJ 
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look in the DMRG measurements to detect the onset of 
the DBL[3,0]. As with the DBL[4,2], the specific loca- 
tions of these features evolve as we vary K/ J and J^/ J 
within the DBL[3, 0] region. As usual, within our the- 
ory, this evolution is a result of increasing ring exchange 
inducing an increase in hopping anisotropy between the 
di and ^2 partons (in DEL [3, 0], the d2 configuration is a 
fixed fully filled ky = band, so only the di configuration 
varies), which in turn results in different "2A;f" wave vec- 
tors as detected by a measurement of _Df,(q). We demon- 
strate this phenomenon clearly in Fig. |8j where DMRG 
data for Di,{qx,qy = ±27r/3) is plotted for various values 
of ring-exchange coupling. 

Another test considers the entanglement entropy. For 
the DEL [3,0], we expect only two gapless modes and 
hence it is computationally feasible to extract the central 
charge from the scaling form of the entanglement entropy 
for comparison. Figure 4 of Ref. 19 shows the fit value for 
the central charge, c, for a range of values of K/ J for fixed 
Jj^/J = 1 on a 24 X 3 system. The jump from c ~ 1 to 
c ~ 2 at roughly the same value of K/J at which the DEL 
phase onsets in the phase diagram is striking. The inset 
of Fig. 4 in Rcf. 19 shows the actual entanglement entropy 



and fitted curves for both the superfluid and DEL phases 
as a function of subsystem size X. In this paper, we 
discuss finite-size effects on entanglement entropy scaling 
below in Sec. IV E (see Fig. 
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For larger values of J_l/J, the system goes into a con- 
ventional rung Mott insulator phase, so called due to the 
decoupling of the rungs; in the caricature of this phase, 
each rung contains exactly one boson which is in the zero 
y-momentum state. These points are marked with trian- 
gles in Fig. [T] Our determinantal wave function actually 
does an excellent job modeling this phase by filling only 
a single band for each parton. This phase has a gap to 
all excitations. We note that very recent work on the 
two-leg XY model [Eq. (11) with K ^ 0] has surpris- 
ingly revealed that at half-filling the superfluid phase is 
weakly unstab le to such a rung Mott insulator for any 
finite J.L.'S^ 

There are two regions, one large above the DEL [3,0] 
and one small at intermediate K/ J and very small J±/J, 
filled with white in the phase diagram. The larger re- 
gion has DMRG points marked as squares with ques- 
tion marks. The VMC data suggest that the DEL[3,0] 
wins energetically in this region, but the DMRG data 
show that the identifying features fail to persist; i.e., 
the singular peaks in the density- density structure fac- 
tor smoothen and the central charge falls below c = 2. 
An analysis of the spectral gap at various values of 
using ED and DMRG reveals that the system is either 
gapless or has a very small gap in this parameter regime, 
as we demonstrate explicitly in the Supplemental Mate- 
rial of Ref.llH Also, we have not been able to identify any 
obvious ordering in the ground state within this region. 

Finally, there is no phase separation for this system. 
The particle density of p = 1/3 is highly stable due to 
the fact that it is commensurate with the number of legs. 



B. The DBL[3,0] phase 
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FIG. 8: (color online). Density-density structure factor eval- 
uated at Qy = ±27r/3, i.e., Db{qx,<ly = ±27r/3), for various 
K/J and fixed Jx/J = 1 on a 30 x 3 system with Nb — 30 
(p = 1/3). The qualitative change in this function is striking 
as one crosses the first-order phase transition from the super- 
fiuid to the DBL[3,0] at around K/J ~ 2.2. In particular, 
Db{qx,qy ~ ±27r/3) is a featureless function in the super- 
fluid phase, but has singular structure in the DBL[3,0] at 
wave vectors originating from the "2^^" wave vectors within 
the di parton band filling configuration. Furthermore, the 
evolution of the singular features within the DBL[3,0] phase 
is fully consistent with the three di parton bands becoming 
more equally occupied as K is increased. For very large K, 
e.g., K/J = 20, the system is in a phase with an approxi- 
mately conserved number of bosons in each rung and chain, 
which is well represented by three equally filled di bands. 
These calculations were done with DMRG. 



We focus now in greater detail on the DEL[3, 0] phase, 
which is an incompressible, gapless Mott insulator (GMI) 
that shows strong density-density correlations at finite 
transverse momenta and at incommensurate longitudi- 
nal wave vectors. We choose a characteristic point, 
K/J = 2.7 and J_l/J = 1, within the GMI region of 
the phase diagram for careful comparison of DMRG and 
VMC results on a 24 x 3 system. For the VMC wave 
function, we again use antiperiodic boundary conditions 
on the partons in the longitudinal (x) direction, but this 
time take periodic boundary conditions in the transverse 
(y) direction such that ky — 0, ±27r/3. The di and d2 
occupations for the VMC state at this point are shown 
visually in the bottom panel of Fig. [9j the ky = band 
for the d2 partons is fully filled with a non-zero gap to 
the other two bands, while the di partons partially oc- 
cupy all three bands. The ky — band has 12 particles 
and the ky = ±27r/3 bands (which are always degener- 
ate due to the inversion symmetry of the lattice) have 6 
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particles each. Therefore, there are no Fermi wave vec- 
tors for d2 (i.e., the band is dispersionless) and we have 
k^°l = 1277/24 = 7r/2 and k^pj^''^'^'' = 67r/24 = 7r/4. 

As mentioned earher, the boson momentum distribu- 
tion function (top panel of Fig. [o]) is featureless, the VMC 
results not showing any structure at all. This is because 
the condition of one particle per rung is exact in the de- 
terminantal wave function whereas, in the DMRG mea- 
surements, the number of particles per rung does fluctu- 
ate slightly. The density-density structure factor (mid- 
dle panel of Fig. [9|, on the other hand, shows singular 

±27r/3) = ±(37r/4,±27r/3) 



features at ±(fc^^^ -I- k 



J±2rr/3) 
Fl 



and ±(2/^^="^""'^^T27^/3) = ±{n/2,T2n/3). 



two potentially enhanced wave vectors with n,y 



±{2kf^ 



The other 

ky = 0, 

p,,0) = (77,0) and ±(4f"/^^ +4f-/3),o) = 
±(77/2, 0), are suppressed by the d2 partons. In the VMC 
data, this is clear due to the fact that the single boson 
per rung condition yields Db{qx, = 0) = exactly. The 
DMRG results show small fluctuations away from zero in 
this quantity with dependence near — 0. This is 
evidence of the phase's incompressibility and rules out 
the possibility of it being some sort of superfluid, either 
conventional or paired. Also, a finite-size study of one- 
and two-boson gapji^l strongly suggests that the realized 
phase is indeed insulating, as does a direct measure of the 
compressibilitjf^ (k — )■ 0) obtained by a scaling analysis 
of bipartite number fluctuations (data not shown) . 

In this paper, we go beyond the results presented in 
Ref. [19iin three ways. First, we investigate flnite-size ef- 
fects in the DEL [3, 0] phase with respect to scaling of the 
entanglement entropy and robustness of the singular fea- 
tures in the density-density structure factor; these results 
are highlighted in Figs. 10 and 11 Second, we make de- 



tailed comparisons of ED and VMC results, focusing on 
the ground state energy and momenta as obtained from 
both methods; these results are highlighted in Figs. [12] 
and 13 Finally, in the limit J = 0, we map our three-leg 



J-J±-K model at p = 1/3 filling to an equivalent spin- 
1 model and suggest a potential connection between our 
c = 2 gapless DBL[3, 0] phase and the c = 2 gapless phase 
known to exist in the bilinear-biquadratic spin-1 chain. ^® 
In Fig. [lOj we plot for various the bipartite 
entanglement entropy, S, versus log(c?), where d = 
{L^/tt) sm{TTX/Lx) is the conformal length and X is the 
subsystem length. According to Eq. (21), at a given 



point in the phase diagram we should observe data col- 
lapse onto a linear function with universal slope given 
by c/3. At the point K/J = J±/J = 1 deep within 
the superfluid phase, the collapse is excellent for all sys- 
tem sizes Lx = 12,18,24,30 onto a curve S{X,Lx) — 
(1.00/3) log(d) -I- 1.12, which is strongly suggestive of 
c — 1 gapless mode as we should expect in the super- 
fluid phase. On the other hand, at the point K/J — 2.7 
and J±/J — 1, i.e., the representative DBL[3,0] point 
presented above and in Ref. IHl the data collapse is less 
impressive, but is still very suggestive of the c — 2 gap- 
less modes predicted by our gauge theory description of 
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FIG. 9: (color online). Boson momentum distribution func- 
tion (top panel) and density-density structure factor (mid- 
dle panel) for a characteristic DBL[3,0] point: J±/J = 1, 
K/J = 2.7. The system size is 24 x 3 with Nt = 24 (p = 1/3). 
DMRG data are joined with solid curves while the VMG data 
are joined with dashed curves, (bottom panel) Schematic rep- 
resentation of the energy-optimized VMC state for this point, 
shown in an analogous way to the DEL [4, 2] state in the bot- 
tom panels of Figs. [5] and [6] Here, the di configuration re- 
sembles a traditional gapless parton Fermi sea, while the d2 
configuration is gapped due to the commensurate p = 1/3 
density. This figure has been reproduced from Ref. |19| 
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FIG. 10: (color online). Entanglement entropy scaling for 
different system sizes on the three-leg system at p = 1/3 in the 
superfluid {K/J = 1) and DBL[3,0] {K/J = 2.7) phases, as 
obtained by DMRG with fully periodic boundary conditions; 
J±/J = 1 in both cases. The dashed lines are fits to the 
scaling form, Eq. | |21[ ). For the superfluid point, all the data 
collapse extremely well onto a curve with c = 1.00 and A = 
1.12, which is the lower dashed line. The dashed line for the 
DEL [3, 0] data is that used in Fig. 4 of Ref.[T5]in which 24 x 3 
data was used and the four smallest/largest X values were 
discarded in the fit; the obtained fit parameters are c = 1.96 
and A = 1.41. The collapse of the ~ 24 and Lx ~ 30 data 
is reasonable, although there is a very weak downward shift 
of the slope in the latter case for the largest X. Shell-filling 
effects or other corrections to the scaling form (which we have 
not considered) are likely the cause of the less impressive data 
collapse in the DBL[3, 0] as compared to that observed in the 
superfluid. 



DBL[3,0] (see Appendix A 2). In fact, we are still likely 
observing some finite-size "shell-filling" effects up to sys- 
tems of size 30 X 3 which makes finite-size scaling of any 
quantity difficult, including the entanglement entropy. 
We note that such shell-filling effects are consistent with 
our parton theory of the DBL[3, 0] phase, and such effects 
actually support our identification of the phase found in 
our model as a DBL[3, 0]. Similar shell-filling effects were 
also seen in Ref. [T3] within the two-leg spin Bose metal 
phase when investigating entanglement entropy scaling 
in that system. 

Obtaining the highly converged entanglement entropy 
data necessary for convincing determination of c becomes 
exceedingly difficult for larger system sizes, especially 
when using fully periodic boundary conditions (PBC) 
as we have employed to this point. We have thus also 
considered OBC in the x direction (cylindrical boundary 
conditions overall), in which case it is much less com- 
putationally costly to obtain very accurate DMRG data 
of any observable on much larger systems. For example, 
at the DBL[3,0] point considered in Fig. 10 on a 24 x 3 
system, retaining D = 4000 {D = ISOOjstates in the 




FIG. 11: (color online). Density-density structure factor at 
qy = ±27r/3 in the DBL[3,0] phase {J±/J = I, K/J = 2.7) 
for various system sizes, L^, x 3, using DMRG with cylindri- 
cal boundary conditions. The main power law singularity (see 
inset) is robust up to systems of size 72 x 3. This is consis- 
tent with the long-distance properties of the gapless DBL[3, 0] 
phase, although we have not pursued a full characterization of 
the long-distance power law behavior. There is some small ir- 
regular behavior of the main singularity on the scale of 2-k/Lx , 
but this is to be expected with cylindrical boundary condi- 
tions and shell-filling effects present in our DBL[3,0] theory. 



DMRG gives a density matrix truncation error of 10 ^ 
(10^1°) with P BC (O BC). As we also experienced with 
DBL[4, 2] in Sec.ImB) when using OBC within DBL[3, 0] 



we observe subleading but apparent oscillatory correc- 
tions to the usual CFT scaling form of the entanglement 
entropy [Eq. (21)], which are to be expected given the 
finite Bose wave vectors present in the density-density 
structure factor. However, for DBL[3, 0] the amplitude of 
these oscillations is significantly weaker than in the case 
of DBL[4, 2] and practically nonexistent for the largest 
subsystem sizes X ~ Lx/2. Then, ignoring the oscil- 
latory piece, and fitting the data to the usual form for 
OBC, S{X, L,) = (c/6)log[(2L,/^) sin {7rX/L,)]+A', we 
find a robust c ~ 2 up to systems of size 72 x 3, with no 
evidence of a downward shifting c which could signify the 
eventual opening of a small gap. 

To provide further evidence for stability of the 



DBL[3,0] phase, we plot in Fig. 11 the density-density 
structure factor at qy — ±27r/3 on several large system 
sizes using cylindrical boundary conditions. As explained 
above and in Ref. .19, this measurement detects the al- 
lowed "2fci?" wave vectors within the di parton band- 
filling configuration and should contain power law sin- 
gularities at various incommensurate wave vectors. In- 
deed, up to systems of size 72 x 3, the dominant singu- 
larity (see inset of Fig. 11) does not smoothen as Lx is 
increased, indicating that the real-space measurements 
look like a power law on these system sizes. Although 
there is some irregular behavior of the singularity on the 
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order of 2-k/Lx momenta, this is likely due to shell- filling 
effects and our use of an open boundary condition in the 
X direction. We caution that it is still possible that the 
DBL[3,0] is weakly unstable in our model on very long 
length scales beyond those considered here, although we 
have not observed evidence for such a trend. As with 
DBL[4, 2], the long-distance properties considered above 
can perhaps be investigate d furt her using entanglement 
renormalization techniques ,^^^1^] which we will leave for 
future work. 

In Appendix I A 2[ we present a bosonized gauge theory 
analysis of DBL[3, 0] and argue that it is, at the very 
least, a potentially stable quantum phase. However, as 
we explain in Appendix |A 2[ there are instabilities out 
of DEL [3, 0] that lead to a fully gapped theory but that 
are difficult to characterize because one cannot construct 
any local observables that correspondingly obtain a fi- 
nite expectation value. The rung Mott insulator is one 
such featureless phase, but it is not the only possibility 
since the sign of the Cooper channel interaction in ques- 
tion may lead to different states. It is conceivable that 
the uncharacterized (white) region above the GMI in the 
phase diagram (see Fig. [?]) is in a fully gapped phase 
(topologically) distinct from the conventional rung Mott 
phase; however, we have been unsuccessful in identifying 
it as such. 

We have also performed a detailed comparison of VMC 
and ED results for the ground state momentum and en- 
ergy using small system sizes. First, we would like to em- 
phasize that ED alone on small clusters can ultimately 
never have the final say regarding the stability of gapless 
Bose metal-like states, e.g., the d-wave Bose liquid or spin 
Bose metal, in a given model. This is clearly due to an 
insufficient resolution of the Brillioun zone for detecting 
a critical, singular surface, but also because the presence 
of incommensurate wave vectors makes finite-size scaling 
difficult to interpret, especially on the system sizes ac- 
cessible to ED. For these small systems, we do expect 
finite-size effects to be significant in our DBL[3, 0] phase, 
and thus the inability to symmetrically fill bands in the 
VMC wave functions to be more prevalent. If these vari- 
ational wave functions truly do capture the main physics 
of the DBL[3, 0] phase, the resulting finite momentum in 
the ground state should be detectible using ED methods. 
In Fig. |12[ we show the ground state momentum com- 
parison for a 10 X 3 system with 10 bosons. Although 
the VMC wave functions do not quantitatively repro- 
duce all transitions between different regions of ground 
state momentum, it is clear that the ED data are re- 
flecting the finite-size effects as manifested in the ground 
state momenta of the nearby GMI states. The particular 
system size that we are presenting actually refiects the 
best agreement over a range of sizes considered (for sys- 
tem sizes Lx X 3, we considered — 4, 5, 6, 8, 9, 10, 12). 
Taking these data sets as a whole, we can only say that 
the ED measurements are sensitive to the finite ground 
state momenta and show some agreement with VMC re- 
sults, while the predictive power of the VMC data for 
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FIG. 12: (color online). Ground state momentum diagram 
for a 10 X 3 system with Nb = 10 (p = 1/3). The colored 
regions denote the ground state momentum of the variational 
wave function that optimizes the trial energy at each point 
in parameter space. The symbols indicate ED data for the 
true ground state momentum of the system. The notation for 
the quantities in parentheses is {Qx,Qy), where Qx and Qy 
are measured in units of 2n/Lx and 2Ti/Ly, respectively. For 
example, (4, 1) ^ (4 x 27r/10, 1 x 27r/3) = (47r/5, 27r/3). The 
agreement between the ED and VMC predictions of ground 
state momenta at this system size is rather remarkable. 



the boundaries between various momentum eigenstates 
in the finite system size corresponding to different band 
fillings is not always compelling. Since we expect our 
DBL predictions to be more accurate as finite-size effects 
are diminished, these results are perhaps not surprising. 
Finally, we note that a recent, analogous consideration of 
momentum quantum numbers in ED ground states tan- 
talizingly points to the possibility of a spin Bose metal 
state existing near the Mott transition in the half-filled 
triangular lattice Hubbard model,^ which serves as a 
nice example of using ED to at least suggest the poten- 
tial stability of a gapless Bose metal-like phase. 

We also performed a study of the J = case, where 
the condition of one particle per rung is strictly enforced 
by the model itself, using VMC and ED. In this ex- 
treme limit, there is no superfluid phase and we expect a 
(nearly) direct transition from the rung Mott insulating 
phase to the DBL[3, 0] phase as we increase the single 
control parameter K/ J±. However, as discussed above, 
DMRG on larger systems indicates that there is likely 
an intermediate phase between the rung Mott insulator 
and DBL [3, 0],'^^! a phase which we have been yet unable 
to characterize fully; this phase also appears at J = 
in the DMRG at larger systems. In Fig. [13) we show a 
comparison of the VMC and ED ground state momenta 
as we vary this parameter for a 10 x 3 system with 10 
bosons. The agreement is excellent over the entire range 
of values studied suggesting that, in this limit, the VMC 
wave functions capture the main physics exceptionally 
well, even though the DMRG and VMC do not agree as 
favorably in the transition region on larger system sizes. 

Finally, in this limit of J = 0, we have also considered 
mapping our three-leg model at p = 1/3 to an equivalent 
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FIG. 13: (color online). Three quantities are plotted on these 
axes for a 10 x 3 system with A^i, = 10 (p = 1/3) with J = 
where K/J± is the only dimensionless parameter. We plot 
the independent variable on a log scale to better interpolate 
between the two limits. Small K/ J± corresponds to the rung 
Mott phase where the rungs of the ladder effectively decouple 
and there exists a gap to all excitations; large K/ J± should be 
the gapless Mott insulator phase, or DBL[3, 0]. The solid and 
dashed curves are the VMC results, while the points are ED 
results; the vertical dashed line is the rung Mott-GMI phase 
boundary as determined by VMC. The blue curve with blue 
circles is the a;-component of the ground state momentum 
measured in units of 27r/10; the red curve with red squares is 
the y-component of the ground state momentum measured in 
units of 2n/3; finally, the black curve with black triangles is 
the ground state energy per site. Remarkably, the variational 
wave functions capture the energetics exceedingly well for the 
rung Mott and the DBL[3,0] phases, much more so than for 
the DBL[4,2] and other conducting DBL[n, m] states. 



model of interacting spin-1 degrees of freedom living on 
a ID chain. Working as always with periodic boundary 
conditions in the y direction, the state of each rung can 
be labeled by the y momentum of the single boson on that 
rung: let |0)i, denote the qy = 0, +27r/3, — 27r/3 

states of rung i, respectively. Our three- leg J±-K model 
at p = 1/3 can then, up to a constant, be written in 
terms of spin-1 operators as follows: 

H = if ^ + Ag,,,+i] + 3Ji (22) 



where we have defined 

P,j{e) = V2 [cos6'(S, • Sj) + sme{S, ■ Sj)"^] 



Q=l,2,3 



with 
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and introduced the parameters A and 9 to connect to 
previous studies of similar models; our J±-K model cor- 
responds to 6* = 7r/4 and A = 1. The operators Sf'^'^ are 
the usual spin-1 operators for the pseudospin on rung i. 

Up to a constant, Pi.j{9 — 7r/4) is the SU(3) symmet- 
ric spin-exchange operator for spin-1. Hence, Eq. (22) 
with 6 = 7r/4, A = 0, and J_l = is equivalent to the 
SU(3) symmetric Heisenberg model, i.e., the spin-1 Lai- 
Sutherland model,!^^^ ^ wh ich is known to be critical with 
c = 2 gapless modes. fS^l^i^tSlI There is also known to be 
an extended critical phase in this model, still at A = 
and Jl — 0, but away from the SU(3) point in the pa- 
rameter range 6 G [7r/4, ^'^ Interestingly, this 
critical phase has soft modeJ^ at — 0,±27r/3 with 
dominant spin quadrupolar correlations at wave vector 
27r/3,'^ which is precisely the "2kp" wave vector in the 
J = J± = limit of our DBL[3,0] theory (i.e., three 
equally filled di bands) at which we should expect singu- 
larities in, e.g., Db{qx,qy — ±27r/3). It thus seems plau- 
sible that the well-known extended critical phase in the 
bilinear-biquadratic spin-1 chairP^l ig connected to the 
decoupled chains limit of DBL[3, 0], although we can't 
be sure without a direct study of Eq. (22) at variable 
A. Finally, the SU(3) symmetric Lai-Sutherland model 
{X ^ Q, 9 ~ 7r/4) has also been studied in the presence of 
single-site anisotropy,'^ which is precisely our J± term. 
As in the pure SU(3) symmetric case, this model is in- 
tegrable via the Bethe ansatz, using which it was shown 
that the c = 2 gapless SU(3) phase exists at finite J± 
before a transition to a large Jl gapped phase, which 
corresponds in the boson language to our rung Mott insu- 
lator. This is consistent with our finding that DBL[3,0] 
is stable for finite values of J±. 

All that being said, our actual J±-K model maps to 
Eq. ([22| at finite A = 1 and fixed 9 = n/A. While the 
operator Qij is rather unnatural for a spin system and 
we have not pursued a direct study of Eq. ( 22 ) with 
variable A, it seems reasonable that the known c — 2 
gapless phas^^ present in the extended parameter space 
discussed above could be stable at A = 1 and thus adia- 
batically connected to our DBL[3,0]. However, we note 
that, unlike the gapless phase in the bilinear-biquadratic 
spin-1 model, DBL[3, 0] at finite J± generally exhibits 
correlations with incommensurate wave vectors. 

Rather remarkably, we have been able to make a (po- 
tential) connection to an exotic spin-1 systenP^ coming 
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from a considerably distant starting point, i.e., itiner- 
ant bosons hopping with ring exchange on a three-leg 
ladder. The above considerations elucidate that it is ul- 
timately the y coordinate of the bosons responsible for 
the gaplessness of DBL[3,G]. In fact, the situation of 
entering DBL[3,0] out of the superfluid by increasing 
K is somewhat reminiscent of interactions inducing a 
metal-insulator transition in a half-filled Hubbard chain, 
wherein the charge degrees of freedom get immobilized 
while the spin degrees of freedom remain gapless. In the 
case of the superfluid-DBL[3, 0] transition, the x coordi- 
nate of the bosons becomes localized, while, very loosely 
speaking, the y coordinate remains free — a more rigor- 
ous interpretation of the phase in terms of bosons only 
is unavailable, hence our fermionic parton description. 
Ultimately, we find it quite profound that the number 
of gapless modes actually increases by one upon enter- 
ing the insulator in the latter case, whereas it always 
decreases in the well-known former case. 



C. Discussion 

Again, the qualitative agreement of the VMC and 
DMRG results for the DBL[3, 0] is excellent. The quanti- 
tative agreement, in fact, is better than in the case of the 
DBL[4, 2]. Furthermore, our confirmation of the central 
charge of c ~ 2 for the DBL[3, 0] and the compelling re- 
sults of the VMC /ED comparison for small system sizes 
serve to reinforce our proposed realization of this phase in 
our model. We also carefully considered finite-size effects 
using DMRG, and, although we can never rule out even- 
tual small gaps appearing on very long length scales,'^ 
we believe the data presented strongly supports identifi- 
cation of the remarkable phase found in our ring model 
as an exotic gapless Mott insulator. 

Subsequent to our discovery and characterization of 
the DBL[3, 0] phase, we revisited the two-leg ladder and 
the possibility of realizing the analogous phase at half- 
filhng {p = 1/2), i.e., a DBL[2, 0] phase. Unfortunately, 
only the most mundane incarnation of this phase is re- 
alized for Jj_ = and sufficiently large K/J. For any 
finite J_L, the DBL[2,0] phase is unstable toward either 
a (tt, tt) charge density- wave (CDW) (for moderate K/J 
and small J±/J) or a state with staggered currents on 
the rungs (for larger K/J). We discuss the situation on 
two legs in greater detail in Appendix [C| 

We also considered the possibility of finding a 
DBL[4,0] phase on the four-leg ladder at either p = 1/4 
or p = 1/2. At filling p = 1/4, DMRG indicates that 
the system phase separates into p — Q and p = 1/2 
filled regions immediately outside of the superfluid, so 
no GMI-type phase is observed. On the other hand, at 
p = 1/2, we observe (tt, tt) CDW ordering even in the 
DBL[4,0] VMC wave functions, thus giving little hope 
for realizing such a phase in the DMRG. However, the 
DMRG itself does not show robust long-range CDW or- 
der even at very large K and reasonably large systems 



(e.g., at K/J = 64 up to systems of size 48 x 4). In- 
stead of a clear Bragg peak, density-density correlations 
[see Eq. ( 16 1] appear to decay as an oscillatory power law 
with wave vector (tt, tt). This same behavior is also ob- 
served in our four-leg half-filled model with unfrustrated 
ring exchange {K < 0) on the same system sizes (e.g., at 
K/J = —64 up to systems of size 48 x 4). Recall that for 
the if-only model (J = 0) the sign of the ring-exchange 
term can be changed by dividing the square lattice into 
four sublattices and performing the canonical transfor- 
mation 6 — > — 6 on one sublatticc, hence the sign of K is 
irrelevant in this limit. We believe the lack of long-range 
CDW order in our half-filled four-leg system, for both 
signs of K, is possibly related to the unusual size depen- 
dence of the (tt, tt) CDW order parameter observed in the 
original studies of the 2D u nfrustrated J-K model with 
quantum Monte Carlo.'^^'^^l For example, at K/ J = —64 
in 2D, the CDW order parameter shows distinct non- 
monontonic behavior with 1/L, and extrapolations per- 
formed using only system sizes up to 8 x 8 would lead to 
the ultimately incorrect concl usion that the system lacks 
long-range (tt, tt) CDW order.l^SIii! We conjecture that a 
similar effect is occurring in our four-leg system due to 
the small transverse size of the ladder. 

This serves as a reminder that the four-leg system is in 
reality still not 2D, and that we should always be cautious 
about extrapolating results from N — 2,3,4, .. . legs to 
full 2D. Exactly what these results at p = 1/2 on four 
legs, which fail to connect particularly well to the known 
2D unfrustrated results, mean for the eventual stability 
of metalhc DBL phases in fuU 2D at p < 1/2 (e.g., 2D 
extensions of DBL[4, 2] from Sec. HI I is presently still un- 
clear and deserving of future investigation, although the 
final answer may not emerge until there exists a scalable 
numerical method suited for L log L scaling of entangle- 
ment entropy in 2D. However, the DBL [4, 2] at density 
p = 5/12 < 1/2 on four legs seems very robust in our 
model, and the behavior of the model at p = 1/2 for 
K ^ J is a special case not necessarily related to the 
stability of the DBL, which is expected to exist only for 
p < 1/2 and intermediate K/J. 



V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In Ref. [IHl the authors presented compelling evidence 
for the existence of a quasi- ID descendent of the ex- 
otic, strongly correlated, two-dimensional d-wave Bose 
metal phase. The so-called DBL[2, 1] exhibited the dis- 
tinctive singular fingerprints of the parent 2D phase and 
possessed the two gapless ID modes (same as the num- 
ber of legs) predicted by the bosonized gauge theory. In 
this paper, we have taken a significant step toward the 
2D limit by finding and characterizing a stable DBL [4, 2] 
phase for our ring-exchange model, which has five gap- 
less modes, one more than the number of legs. In par- 
ticular, both of the slave fermions, di and c?2, occupy 
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more than a single momentum band and thus manifest 
the d-wave sign structure of the wave function in a highly 
nontrivial manner, demonstrating extensibility to two di- 
mensions. It may yet be possible to consider a system of 
six legs {N = 6), using the results obtained thus far for 
= 2,4 to guide the search for perhaps a DBL[6, 2] or 
even DBL[6, 4] phase. The great challenge here is that 
the spatial entanglement becomes distressingly large as 
the gaplessness increases and DMRG methods will have 
great difficulty converging meaningful results on systems 
of reasonable length. Also, as we discuss further in Ap- 
pendix |B] it may be possible to stabilize DBL phases on 
odd leg ladders with p > 1/N by modifying the trans- 
verse boundary conditions; this would potentially allow 
us to realize DBL phases on three- and five-leg ladders 
that are also extensible to 2D. [In the case of the three- 
leg ladder with periodic transverse boundary conditions, 
our numerics indicates that such Bose liquid phases are 
not stabilized, and instead an exotic ("bond chiral") su- 
perfiuid phase'^'^ is likely realized by our ring model (see 
Appendix B 2 for details).] 

Our tried and true method of approaching 2D via a 
series of quasi- ID ladders has remarkably yielded great 
success in demonstrating the DBL as a potentially sta- 
ble, gapless, and strongly interacting state of hard-core 
bosons for which no perturbative understanding cur- 
rently exists. The quasi-lD nature of the systems we 
have studied has allowed us to access a very challenging 
realm of condensed matter physics both with exact nu- 
merical techniques (DMRG/ED) and with a robust gauge 
theory that, in principle, allows for the exact calculation 
of all anomalous power law exponents. It is worth not- 
ing that it is precisely the existence of singular surfaces 
in momentum space and the lack of a quasiparticle de- 
scription in the 2D phase that simultaneously render the 
numerical and analytical analysis extremely formidable 
while paving the way for our quasi- ID methods to be es- 
pecially effective in gaining a physical intuition for these 
strongly correlated phases. It is only by virtue of the 
"Bose surfaces" that we can reliably and meaningfully 
detect the corresponding "Bose points" on the ladders. 

However, we still caution that the ladder approach is 
not perfect, and studies at N = 2,3,4,... legs hardly 
represent crossing over to full 2D. Indeed our results on 
ladders of this size are somewhat irregular. As discussed 
further in Appendix |B 2[ we do see a strong even-odd 
effect (with the number of legs N) in the stability of 
compressible Bose-metal phases; the transverse bound- 
ary condition is a substantial modification to the model 
on these few-leg ladders, and we conjecture that altering 
it could change the story considerably. We note that in 
some cases it is possible to truly solve frustrated 2D mod- 
els with ingenious and systematic use of the DMRG,1S2MI 
but unfortunately, those techniques are not practical for 
use with our Bose-metal-type phases in our ring model. 
Nonetheless, the local physics is expected to become more 
2D-like as N increases, and thus we believe our four-leg 



tual stability of the DBL in 2D. 

In the course of our search for gapless d-wave Bose 
metals, we came across the highly novel gapless Mott in- 
sulating phase at commensurate density of p = 1/3 on the 
three-leg ladder. The so-called DBL[3, 0] gave rise to os- 
cillations at incommensurate wave vectors in its density- 
density correlator, possessed two gapless modes, and ex- 
hibited the characteristic d-wave correlations common to 
all of the DBL[n, m] phases. Despite being truly quasi- 
lD in nature, this strange insulator is very interesting in 
its own right as it can be described by a parton gauge 
theory in which one parton is gapless while the other is 
gapped (see Appendix A 2). We expect that it may be 



results in Sec. Ill are at the least very suggestive of even- 



possible to find a quasi-2D analog of this phase by study- 
ing an N-layer square lattice system at density p — j /N 
for j — 1,2,... where we believe that strongly coupled 
ring exchanges between the layers may give rise to the 
desired behavior. Additionally, while our study of the 
four-leg ladder did not discover a stable DBL[4, 0] phase, 
it may be the additional frustration of an odd number of 
legs that allowed the DBL[3, 0] to exist stably. Therefore, 
it may be possible to search for an analogous DBL [5,0] 
phase on the five-leg ladder at densities p — 1/5,2/5. 
The smaller density will likely phase separate right out 
of the superfluid phase as was the case for p = 1/4 on 
the four-leg ladder, but the larger one may not. 

Finally, we should here again mention the physical 
problem that inspired the original search for the DBL 
in the first place: itinerant electrons with real-space d- 
wave correlations; that is, a d-wave metal phase. If one 
imagines fractionalizing the electron creation operator 
into a charge-carrying hard-core bosonic piece and a spin- 
carrying fermionic piece (i.e., cj. — fl) and then frus- 
trating the system sufficiently (with an appropriate ring 
term) such that the bosons fail to condense and instead 
realize a DBL phase (i.e., — d\d\), the resulting elec- 
tronic metal should be accessible by all of the methods 
employed herein. Namely, we can place such a system on 
a ladder and study it numerically using DMRG and ana- 
lytically using a similar bosonized gauge theory approach 
as that presented in this work. Such an investigation is 
currently underway with promising results thus far. 
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Appendix A: Gauge Theory Description and 
Solution by Bosonization 

Here we describe the analytical approach for the A^-leg 
ladder where > 3, which is generic; the formulations 
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for iV = 1, 2 are special cases. To accurately capture the 
behavior of the physical bosons within the slave parti- 
cle picture [see Eq. ^] requires a compact U(l) lattice 
gauge theory. First, we consider the gauge sector: we 
denote the integer-valued electric field on links between 
sites X and x + 1 on the same chain of the ladder by 
Ei ;^, where £ — 1, . . . ,N labels the chains, and on links 
between chains i and i + 1 aX the rung located at x by 
Ei,y. We then introduce a 27r-periodic vector potential 
on each link of the lattice, which is denoted the same way 
as the electric field: ae,x for links within chain i and a^ y 
for links between chains £ and £ + 1. The magnetic field, 
which lives on the elementary plaquettes is then given by 
the lattice curl of the vector potential: 

(Vxa^) ^ ai^x{x)+ae^y{x+l)-ae+i^x{x)-ai^y{x), (Al) 

where ai — {ai^x,o-i,y)- The Hamiltonian density for the 
gauge sector is thus 



N 

E 

e=i 



4 



K cos{V X ae) 



(A2) 



where = {Ei,x, Ee^y) and where we identify £ — N + 
1 = 1 due to the transverse periodic boundary conditions. 
We also include here arbitrary coupling parameters for 
the electric and magnetic fields, h and K. Note that 
the contribution due to the magnetic field is encoded in 
a compact form due to the 27r-periodicity of the vector 
potential. 

The purpose of the gauge field in our slave particle 
theory is to bind the di and ^2 partons together at each 
site to realize the hard-core bosons. To this end, we 
assign equal and opposite gauge charges to the partons; 
specifically, charge ei = +1 to the di partons and charge 
62 = —1 to the ^2 partons. The lattice divergence of the 
electric field gives the total free charge: 



(y ■Ei)^Ei^,{x) 



Ei^x{x - 1) + Ee,y{x) - Ee^i^y{x), 

(A3) 



(V • E^) = 4^1 - dld2. 



(A4) 



We therefore wish to enforce the constraint that (V • 
E^) = at each site. This is accomplished by adding a 
constraint term to the action. 

We work with a Euclidean path integral and denote 
the temporal components of the vector potential at each 
site on chain £ as a^^r- The electric field and vector 
potential are canonically conjugate quantum operators: 
[Eg^fj_, ai^v] = iS^j^ij and so the full gauge sector Lagrangian 
is thus 



N 



E? 



V^a^ + — K cos{y X a^)-!- 

+i(V-E<>)af,^], 

(A5) 



where 'S/ T-ai,u = 0'i,v{x, t -\-\) — ai ,^{x, r). The first term 
is the "ipdrQ^^ term where [q,p\ = i. The second and 



third terms are the Hamiltonian. The final term is the 
constraint that will enforce (V • E^) = upon integrating 
out a£_T-. By taking advantage of the 27r-periodicity of the 
resulting action in the spatial components of the vector 
potential, we can choose to integrate out the electric field 
over all real numbers despite the fact that it is strictly 
integer-valued. The resulting Lagrangian is 



N 

E 



e=i 



1 

2h 



2h 



(VrCt^,!^ — ae+l.T + o-e,T)'^ — /'ircos(V x ai) 



(A6) 



where VxO-e,!^ = ai.u{x -f 1, r) — ai^^{x, r). 

We now turn our attention to gauge fixing. In princi- 
ple, we have the freedom to specify a scalar field A^(x, t) 
at every site in the space-time lattice; however, in antici- 
pation of taking the continuum limit in x and r, we only 
have the freedom to choose N functions of x and r, which 
are precisely the functions A^. Let us first consider how 
the a^^y fields transform under a gauge transformation. 



ae^y + Ki+i{x, t) - Ai{x, r). 



(A7) 



Interestingly, the combination ay = (l/A^)[X]fci '^^.y] 
thus inherently gauge invariant. If we imagine our lad- 
der as a cylinder due to the strictly finite extent and 
periodic boundary condition in the transverse direction, 
this combination of gauge fields corresponds to the mag- 
netic field through the cylinder, which we can choose to 
set to zero without expending any gauge freedom. We can 
then choose — 1 constraints on the scalar fields A^ so 
as to fix all the other y-component gauge fields to zero: 
ai,y = for all £. This leaves us with exactly one con- 
straint left to choose. The Lagrangian can be simplified 
at this stage as follows: 



N 

E 

e=i 



1 

2h 



+ ^(a^+i,r - a^.r)^ - Kcos{ai+i^^ - ai^^) 

(A8) 

In formulating the fermion sector of our gauge theory, 
which we will do below, we will consider as independent 
flavors the partons corresponding to different momentum 
bands. These partons couple to the gauge field combi- 
nation a^ = {i/N)[^^^-^ai^^], where /i = x,t, as to 
the usual gauge field. The other linearly independent 
combinations of the gauge fields, ai^^, are all massive in 
Eq. ( A8 ) and can thus safely be integrated out of the ef- 



fective Lagrangian, leaving only Ox and a^. We can now 
use our final gauge constraint to fix a.j. — 0. Thus, the 
final, effective gauge sector Lagrangian is quite simple: 



N 



(V^ar)^ 



(A9) 
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where we have taken the parameter h = 1. 

We now consider the fermion sector, in which we start 
with just the two partons that compose the hard-core 
bosons in the slave-particle picture: da(v) for a = 1,2. 
Then we decompose each such operator into its discrete 
Fourier components in the y direction: 

di{x,y) = -^^e^'=.fe'i)rf('=»)t(^), (AlO) 

ky 

where ky runs over the partially occupied bands for each 
value of a. For the DBL[?i,m], ky varies over n discrete 
values of the transverse momenta for a — 1 and over m 
discrete values for a = 2. From this point on, we take the 
continuum limit in the longitudinal direction and treat x 
as a continuous variable. All operators in what follows 
will be functions of this single variable x and so it will be 
suppressed henceforth. 

Each partially occupied band has orbitals filled on the 
interval [—kp^\k'"p^^] and this defines the Fermi wave 
vectors. We are primarily interested in the low-energy 
physics of the DBL phases and so our next step is to 
split the fermion operators into right- and left-movers 
linearizing about the right and left Fermi points for each 
parton: 



where b, b' are composite indices for a and ky and run 
over all possible values of these. The density operators 
are simply pbp = d\pdbp. The couplings Bb^},i,Ch,b' can 
be thought of as (m + n) x (m -I- n) symmetric matri- 
ces for a DBL[m, n]. Such terms are strictly marginal 
with respect to the DBL fixed point, but will shift and 
renormalize the initial dispersion velocities. 

Other, potentially more worrisome, four-fermion inter- 
actions do exist and will be considered within the context 
of particular DBL[n,r7i] phases below. 

The gauge theory as presented thus far is entirely treat- 
able by conventional bosonization techniques. This con- 
trolled approach allows one, in principle, to compute all 
power law exponents thus capturing all the long-range, 
low-energy physics of the putative DBL[n, to] phases. We 
begin with the bosonization formula 



(AM) 



'■aL ^ 



(All) 



where (j)a^^ and olt^^ are the conjugate phase and phonon 
fields for each band satisfying 

(^^^^\x)/^P{x')\ =i7rS^^,Sk^,,eix-x'), (A15) 

(k ) . 

while Tja " are Klein factors, which preserve the fermion 
anticommutation relations, i.e., they commute with the 
bosonic fields and anticommute among themselves. The 
fermion sector Lagrangian, cast in bosonic fields, be- 
comes 



which gives the full right-left decomposition shown in 
Eq. Throughout this paper, we use the standard 

conventiorP^ that the grouping of a pair of right- and 
left-movers corresponds to a pair of fields related by time 

reversal: {kx,ky) {-k^^-ky). From now on, ky on Cf^^^-,^ = X! ^ [dj'-^yA (dr(t>'i 

the fields is a label: right-movers carry the same y mo- ^.ky ^ 

mentum as the label, while left-movers carry opposite 
y momentum. Each right- and left-moving fermion be- 
haves like a massless Dirac fermion. With the spatial where 
components of the gauge field [ux and Oy) fixed at zero 
by our chosen gauge constraints, the free fermionic action 
is governed by a Lagrangian of the form 



\ ^ V 



a.ky 



(fey) 
Fa 

2tt 



{dxCi 



(ky 



ky) 



(A16) 



(A17) 



where 



form 



,(fey) 



_,(fcy)t 



i.ky,P 



Qj CL' 



-iPv 



(A12) 



are the dispersion velocities corresponding 



to each Fermi point. In a theory free of any four-fermion 
interactions, these velocities are the only input parame- 
ters. We could, in principle, use variational wave func- 
tions, compared against DMRG results, to determine on 
finite-size systems the approximate filling of the bands 
for particular points in our J-J±-K model [sec Eq. ( 11 )]. 
Assuming free fermion band dispersions, we could then 
estimate the dispersion velocities as simply the slopes of 
the band curves at the Fermi points. 

We now consider potential four-fermion interactions, 
starting with density-density interactions: 

^p-int = [Bb^b'PbPPb'P + Cb,b'PbPPb'-p], (A13) 

bM\P 



The first term in the Lagrangian can be broken into 
two, one of which does not involve the gauge field. This 
term is the "ipc^T?" term where we identify q ^ (/) and 
P = {dxS)/TT ^ p; indeed [4>,p] — i- The other piece of 
the first term, which does involve the gauge field, is a 
descendent of the constraint that will end up fixing the 
di and d2 occupation numbers to one another at every 
position. The rest of the Lagrangian is the Hamiltonian. 



Taking now the full action given by £ = C' +C 



form : 



we can integrate out the gauge field Or, which presumes 
the strong-coupling limit, and arrive at 



^' = J2- {^A'"^) [drcblt'^) + M9^c+n^in, (A18) 

a,ky 



where 



1 Tn * ^ 



(A19) 
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with corresponding mass M = M{N, n+m, h), and n and 
m refer to the labels of the DBL[n, m]. The normalization 
factor is chosen in anticipation of performing a unitary 
transformation on the bosonic fields. This overall charge 
mode, Oc, has been rendered massive by the constraint 
placed by the gauge field. We can thus integrate out the 
(0c: ^c) mode effectively pinning 9^ to zero throughout 
the theory. This condition implies that 



2 



Pd2 



(A20) 



where the ky index in each sum varies over only the par- 
tially occupied bands for each flavor of parton. There- 
fore, there will be no free gauge charge in the system (i.e., 
confinement), which is equivalent to saying that either a 
position will be unoccupied or it will have a hard-core 
boson. 

We now have a theory of m + n — 1 gapless ID modes, 
and, after diagonalizing to find the normal modes, the 
effective Lagrangian can be written as follows: 



n+m— 1 



• * TT 



with 



n-^-m— 1 



1^=1 



2tt 



9y 



(A21) 



(A22) 



The quantities are the Luttinger parameters and v^, 
the dispersion velocities for each normal mode. With 
these ingredients, plus the transformation matrices that 
diagonalize the Lagrangian leading to Eq. ( A21 ), one can 



compute any quantity of physical interest (ignoring four- 
fermion interactions). 

Finally, it is worth noting that the density-density in- 
teractions are trivially bosonized using 



PbP 



1 

2^ 



{d.Jb + Pd^cj),,). 



(A23) 



1. DBL[4,2] 

We work on the four- leg ladder (A'' = 4). Motivated 
by the numerical VMC and DMRG results, we conclude 
that for the DBL[4, 2], the di partons fill all four bands 
in such a way that pairs of neighboring bands have the 
same occupation; that is, there are two pairs of coin- 
cident bands. The ^2 partons fill only two bands and 
these two are also coincident. Thus, the most sensible 
choice of transverse boundary conditions for the partons 
is antiperiodic (this still results in the bosons obeying 
periodic boundary conditions) wherein the ±7r/4 bands 
are coincident as are the ±37r/4 bands, consistent with 
symmetry of the lattice. Hence, there are only three 
unique Fermi wave vectors, which we designate as follows: 



«F11 



^(±'r/4) 



k 



F13 



i,(±37r/4) 
"Fl 



, and kF2 



well as the corresponding Fermi velocities vpii 



f.(±'^/4) _ 
F'2 

(±7r/4) 



V 



= „(±^/4) 
^F2 



Fl 



For boson density 



WF13 = Vpi ' and VF2 

/9, the Fermi wave vectors satisfy kpii + ^fis = 27rp and 
kp2 — 2TTp. In the mean-field treatment of the slave par- 
ticle theory, di and c?2 simply fill up energy bands about 
kx — 0; for free fermions with nearest-neighbor hopping, 
these band curves are given by cosine functions and hence 
we can sensibly assume that the band velocities can be 
calculated simply by evaluating the derivatives of these 
curves at the Fermi points. There are only two distinct 
ratios of velocities in general: 



sin (A; 



Fll) 



VFll _ 

VF2 am{kF2) 

VF13 Sm{kF13) 
VF2 



sm{k 



Fll) 



sin(fci?2) 



sin(27rp) ' 
sin(27rp — fc^ii) 
sin(27rp) 



(A24) 
(A25) 



and hence, once one specifies the density p, there is only 
one input to the theory, kFu (ignoring all four-fermion 
interactions). Since we were very successful in finding 
the DBL[4, 2] at p = 5/12, we shall use this density here. 
A thorough analysis of the scaling dimensions of vari- 
ous boson operators, density excitation operators, and 
four-fermion interactions can now be performed with re- 
spect to the DEL fixed point while varying the input 
parameter fc^ii from vrp, which occurs when the four 
bands are equally filled (this limit corresponds to the 
most anisotropic Fermi seas for the partons and is ex- 
pected for large K and/or small J± in our model), to 
2ttp, which occurs when the ±37r/4 bands are no longer 
occupied for the di partons (in this limit, both flavors of 
parton fill the same Fermi sea and the variational wave 
function is simply a determinant squared; the resulting 
phase has no d-wave character at all and is likely unstable 
to a superfluicP^ . 

This analysis confirms that the boson operators ex- 
pected to have enhanced power law behavior (relative to 
the mean field) via the Amperean rule do indeed have 
scaling dimensions less than one for the range of kFii 
given above, which, for the single-boson correlator, re- 
sults in exponents less than the mean-field value of two. 
These operators, listed by the y momentum they carry, 
are as follows: 



% = 

qy = ± 



TT 

±- 

2 

qy = TT 



,(±77/4)1 ,(±7r/4)t . 

,(±7r/4)t ,(T7r/4)t . 
"l_R "2L ' 

,(±37r/4)t ,(±^/4)t . 
"iJ?, "2L ' 

,(±37r/4)t ,(T^/4)t 
"li?, "2L 



(A26) 
(A27) 

(A28) 

(A29) 



We do not list here operators that can be obtained by 
exchanging i? -H> i. 

Next, considering the density excitation operators, the 
bosonization analysis confirms that the scaling dimen- 
sions of the operators expected to have enhanced power 
laws are indeed less than one. These operators, again 
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listed by their y momentum, are as follows: 



Qy 



. ,(±7r/4)t ,(T^/4) . 
„ . ,(±37r/4)t ,(T3^/4) . 



9y = 
9. = ±2 
9. = ±2 

TT 

9. = ±2 

9. = ±2 
=7r 



,(±7r/4)t ,(T^/4) . 
"'2R "■2L ' 

,(±7r/4)t ,(±7r/4) . 
,(±37r/4)t ,(tt/4) . 
j(T37r/4)t ,(T37r/4) . 
,(±7r/4)t ,(±7r/4) . 
,(±37r/4)t ,(±7r/4) 



(A30) 
(A31) 
(A32) 
(A33) 

(A34) 

(A35) 

(A36) 

(A37) 



of such bilinears will have scaling dimensions greater than 
or equal to two. Bosonization analysis indeed reveals that 
these terms are strictly marginal over the range of kpn 
studied. 

While there are certainly several different phases that 
can be realized by the full phase space of quadratic 
bosonized models, we have showed convincingly that with 
only gauge projection and no residual forward scattering 
interactions, our parton model lies entirely within the 
DBL[4, 2] phase. We have also analyzed the scaling di- 
mension of the cosine terms of Eqs. (A38l and (A39) 



Again, we do not list the qx reversed partners of each of 
these terms. 

We have also considered all possible four-fcrmion in- 
teractions allowed by the symmetries of the system (i.e., 
parton number conservation and crystal momentum con- 
servation). There are the density-density interactions of 
Eq. (A13I, which are strictly marginal for any value of 



kpu since their bosonized form [see Eq. (|A23[)] is hn- suits presented in Sec. |IIIB 



while including various density-density interactions [see 
Eq. (A13)] in the Lagrangian. For the interactions con- 
sidered, e.g., attraction between bosons, we again found 
the cosine terms to be at worst marginal; however, we 
caution that this study was not exhaustive, and it is pre- 
sumably possible to add interactions that will render the 
cosines relevant, especially in the decoupled-chains limit 
where Eq. ( |A38[ ) is already marginal in the pure gauge 
theory. All in all, our analysis still strongly suggests that 
DBL[4, 2] is a stable quantum phase over a large range 
of reasonable parameters. These results, combined with 
the encouraging agreement of the DMRG and VMC re- 

suggest that the DBL[4, 2] 



ear in the bosonic fields; therefore, density-density prod- 
ucts only contribute to the quadratic terms in the theory. 
There exist other terms, however, which yield exponen- 
tials of the bosonic fields and therefore interact beyond 
quadratic order. The first class of such terms have the 
following form: 



, _ ,(a7r/4)t ,(/37r/4)t ,(737r/4) ,(A-37r/4) „ 

n-afi^s - difl fliL difl a-iL + ^-C- 



(A38) 



where a, 13, j, 6 — ±1 and where we must have a — /3 — 
87 -|- 3(5 = 8n for n an integer in order to conserve mo- 
mentum. There are thus three terms, , 

and ^-1 ^, along with their qy reversed partners. Ig- 
noring all four-fermion interactions in the Lagrangian, 
we have performed an analysis of the scaling dimension 
of these terms and found that they are strictly irrelevant 
for the entire range of kpn except for at irp (four equally 
filled bands; large K and/or small J±) where they are 
marginal. Finally, one can form valid four-fermion terms 
from combinations of parton transfers on the same side 
of coincident bands: 



l{a,ky.P),{a',k'y,P') — Ba,ky,pBa' ,k'y^P' + H.C. , (A39) 



where 



(A40) 



To conserve momentum, we must have Pky 
with n an integer. Without doing any work, we would 
expect that these terms are not relevant since the Ba,ky,p 
bilinears are precisely the suppressed density excitation 
operators that do not follow the Amperean rule. We 
expect their scaling dimensions to thus be greater than or 
equal to one and therefore that the four-fermion products 



is a stable phase, the qualitative features of which are 
faithfully manifested by our proposed model, Eq. ( [TT] ), 
over a large region of the K/J — Ja^/ J parameter space. 



2. DBL[3,0] 

In this section we work on the three-leg ladder (TV = 3) 
at commensurate density p — 1/3. The qualitative agree- 
ment of the VMC and DMRG data suggests a DBL[3, 0] 
phase where the di partons occupy all three bands, two of 
which are coincident, while the c?2 partons occupy fully a 
single band. This structure suggests that periodic trans- 
verse boundary conditions for the partons are appropri- 
ate such that the ky — -izlir /i bands are degenerate and it 
is the ky = band that the ^2 partons fully occupy. The 
lack of any partially filled ^2 bands implies that the ad- 
dition of a d2 parton is gapped, and the same must then 
be true for the boson in the strong-coupling limit where 
the partons are bound together. This fact reveals itself 
in the featureless boson momentum distribution function 
obtained via DMRG (see Fig. |9] top panel). Addition- 
ally, the single fully filled band enforces a condition of 
exactly one particle per rung on the ladder. This con- 
straint is exact in the VMC wave functions and even in 
the DMRG, the rung occupation number in the putative 
DBL[3, 0] region rarely differs from one. 

The bosonization analysis of this phase is similar in 
many regards to that of the DBL[4, 2] above, but here 
there are no fields corresponding to the ^2 partons since 
there are no d2 Fermi points. This requires a careful 
encoding of the one-particle-per-rung constraint. We de- 
note the Fermi wave vectors in the three di bands as 



«F0 



kpl and kp = k^^^^^^"^ and the corresponding 
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Vpl and Va 



vf^-'/'l At 



Fermi velocities as vq 
p = 1/3, we have the condition, kpQ + 2kp = n. If we 
again assume that the Fermi velocities are given simply 
by the slope of nearest-neighbor hopping band curves at 
the Fermi points, then there is only one distinct ratio of 
velocities given by 



IV _ sin[(7r - fcfo)/2] 
vq sin(fcFo) 



(A41) 



and so, ignoring residual, short-range interactions beyond 
the gauge field, the theory has only one input parame- 
ter, kp(j, which can be varied from 7r/3 where the three 
di bands are equally occupied (i.e., the large K and/or 
small J± limit) to tt where the state becomes a rung 
Mott phase with the ky = Q band filled fully for both 
partons. This setup allows us to compute various scal- 
ing dimensions numerically as we did for the DBL[4, 2] 
above. However, the two-mode harmonic liquid theory 
for DBL[3, 0] is sufficiently simple that it can be treated 
analytically, as we will now explain. 

We do not expect there to be any long-range boson 
correlations since the ^2 dispersion is gapped; further- 
more, the fully filled ky — band for the ^2 partons sup- 
presses any density-density correlations at ky = 0. We 
do, however, expect there to be singular density-density 
correlations at ky = ±27r/3. Also, there are potentially 
relevant four-fermion interactions that must be consid- 
ered in discussing the stability of the DBL[3, 0]. 

The gauge field construction for this phase follows the 
same approach as for the DEL [4, 2] with the gauge con- 



straint requiring that d[di = d2'^2- But since the d2 
fermion fills a band, there is precisely one d2 fermion 
per rung. At longer wavelengths down the ladder, the 
one-dimensional ^2 density does not fluctuate at all. 
Within the gauge theory the effective constraint is thus 



Sfc„,p d,p 



d^e'^'^yyTT = PA = 0. Here, pA 



is the total "gauge charge" density. Note we are drop- 
ping the subscript a on the partons and the bosonic fields 
since this theory only has a = 1. 

The imaginary-time bosonized Lagrangian density is 
given by Eq. (|A18|): 



(A42) 

where we have adopted the shorthand ± o ±27r/3, 9c is 
given by Eq. (|A19|), and Hkin by Eq. (|A17|. 



To proceed, we make a canonical transformation 
1 



V2 



(+) - 



(A43) 
(A44) 



V3 



and identical definitions for the ^^^f) fields. These can 
be inverted as 



0(0) 



1 „ /T„ 



(A46) 
(A47) 
(A48) 



2 " V 6 

We now integrate out the massive field 6c to obtain 



^kii 



E 



27r 



1 

5m 



with g^r = 1, 



5p = 3 



VaVO 



VcrVo + 2{vq + V^)"^ ' 



'vaVoiVg + 2vq) 
Vq + 2Va 



(A49) 



(A50) 



(A51) 



As it stands, gp < 1, but residual short-range forward 
scattering interactions will renormalize both gp and go- 
as well as the two velocities. Quadratic cross terms that 
couple the two sectors, such as dxdpdxOa are precluded 
by y — !■ —y symmetry. Indeed, under this symmetry, 
0(+) O 0(^) so that Op -J> Op is even and 9^ -0„ is 
odd. The cross term dxOpdxOa- is odd, and not invariant 
under this symmetry, thereby being forbidden. Thus, the 
fixed point theory o f the DEL [3, 0] phase is given by the 
Hamiltonian in Eq. ( A49 1 , or in Lagrangian form by 



£'= V -{dxO^)idrq^^) + nM, 

• * TT 



(A52) 



and consists of two gapless bosonic modes. 

There are two allowed four-fermion interactions other 
than density-density type (already absorbed in the gen- 
eral Luttinger parameters gp and g^)', they can be crudely 
viewed as "Cooper channel" interactions: 



Hi — w\ 



■dWt,(+)t4-)4) 



H.c. 



2wi cos(2(/)(+) - 20(")) 
2w\ cos(2\/2(?!)ct) , 



W2 



2w^ 



cos(20('^) -20^^^) 



(0)^ 



H.c. 



4^2 cos(-\/20o) cos{V &4>p 



(A53) 

(A54) 
(A55) 

(A56) 
(A57) 
(A58) 



(A45) By re-expressing these operators in terms of ( 



one can readily deduce their scaling dimensions: Ai 
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2/(7o-, A2 = 3/(2(7p) + l/(2gcr)- Provided these are irrele- 
vant, the DBL[3,0] phase is stable. This requires 5,7 < 1 
for Ai and gp < 3/(4— l/g^) for A2. The b are gauge the- 
ory that gave go- — I and gp in Eq. ( A50 1 automatically 
satisfies these requirements. 

In the DBL[3, 0] phase the boson operator b = did2 
will be short-ranged because the 6,2 fermion has a gap. 
However, power law behavior is expected in the density- 
density correlator. To examine this we consider var- 
ious density bilinears that contribute to b^b ~ d\di. 
We consider first the dominant "2A;f" contributions at 
ky — ±27r/3, listed by their momenta: 



(2fci.,27r/3) : ^ ^ -^^"^ 



which have scaling dimensions 
(2fc^,2V3):| + |; 



(A59) 
(A60) 



(A61) 



(A:Fo + fcf,27r/3) 



1 /^5p 



3 
9p 



9a 



(A62) 



[Our writing here and below is somewhat schematic, e.g., 
we do not spell out all complex pre-factors and Klein fac- 
tors as they do not affect the power-laws. For the {kpQ + 
kp, 2Tr/3), there is another contribution at the same 
momentum, and a more precise expression would read 

d(^+)t40) ^ ^(0)t^(+) ^ _^^(+)^(0)g-.(.(«)+fl<+)) 3i^(^(0) _ 
^(+)).] 



There are also "2fcf " contributions at A:„ = 



(2fc^o,0):drdf 



= e 



y 

i2^2/3ep . 

i(6/(-)+e(+)) 



(A63) 
(A64) 
(A65) 



with scaling dimensions 2gp/?t and gp/(S-\-\/{2g„), respec- 
tively (the second equalities follow from setting 9c = Q 
consistent with the pinning condition). It would appear 
that this leads to a power law decay of density-density 
correlations at ky — 0, but due to the full band for the d2 
fermion we know that this cannot be correct. However, 
there exist four-fermion operators that will contribute to 
the density and that have the same momenta and scaling 
dimensions as the bilinears above: 



(A66) 



{0)t^(0)^(+)t^(7) _ ^i(0<-)-0<+')^j(2e<«)+9< + )+9<-') 



e 



(A67) 



The presence of these operators follows due to the 
existence of a six-fermion umklapp term wherein a 
fermion in each of the three bands is backscattered, 



a more microscopic implementation of the gauge con- 
straint wherein the di fermion is strictly enslaved to the 
c?2 fermion, one would anticipate that these two opera- 
tors exactly cancel one another when contributing to the 
microscopic density. 

Having considered the potential destabilizing interac- 
tions and determined that they are irrelevant when short- 
ranged density-density interactions are ignored, we can 
conclude that the DBL[3,0] phase is potentially stable. 
In light of the DMRG/VMC correspondence shown in 
Sec. IV B[ it is evident that this phase is likely the ground 
state of our model, Eq. ( [Tl] ), over a large region of the 
phase diagram. 

Let us briefly consider what happens when both inter- 
actions in Eqs. (A53l and (A56l become relevant thus 
pinning the field s (j)p and <!)„. For simplicity, let us as- 
sume that Eqs. (A53) and (A56) represent interactions 
when they are already 0(1) after some initial flows, and 
we now need to minimize Hi -\- 7^2 semiclassically. The 
resulting states depend on the signs of the couplings wi 
and W2, and we consider different cases in turn. 

When wi < 0, we can simultaneously minimize Hi 
and ^2. For either sign of W2-, there is a unique physi- 
cal state (analyzed as in Ref. Uni Sec. IV.E.l.): = 
= 0(0) -I- [1 -I- sign(w2)]7r/4. We did not find any 
local physical observable that would obtain an expecta- 
tion value for such pinning of phases for either sign of 
W2, and hence these states are good candidates for the 
featureless rung Mott insulator phase. At this point, we 
do not know whether the states obtained for W2 > 
or W2 < are qualitatively distinct and cannot be con- 
nected by any path; if they are distinct, then we do not 
know how to interpret the featureless state other than 
rung Mott insulator. It could also be that the two cases 
are only quantitatively distinct and connect to the same 
rung Mott insulator picture. 

When wi > and sufficiently large, wi > |?i)2|/2, the 
energy is minimized at (^^"'"•* = ± acos(— W2/2wi), 
(^(0) _ (0(+) -I- 0(~))/2. For any W2, there are two phys- 
ically distinct solutions, and now we have an "order pa- 
rameter" sin(2(/)("'") — 2(/)(~)) that obtains an expectation 
value and takes opposite signs for the two physically dis- 
tinct solutions. The order parameter respects lattice 
translation and inversion symmetries, but changes sign 
under both time reversal and mirror symmetry. A pos- 
sible state with similar symmetries is a chiral state with 
spontaneously generated current circulation, which can 
be realized, e.g., by introducing uniform flux in the di 
hopping problem. An operator with similar transfor- 
mation properties as the current circulation is x(i") ~ 
J(l (r) — Jj^{r + x), whose correlators can be reduced to 
current-current correlators, (x(i")x(r')) = C^'*'(r, r') -|- 
Cf'* {r + x,r' + x)- C*'^ (r , r' + x) - Cf'^ {r + x,r'). Since 
we have not observed long-range order in current-current 
correlations in any phase at p = 1/3, we conclude that 
this phase is not realized in the ring model. Note that 
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FIG. 14; (color online). Phase diagram of the three- leg sys- 
tem at boson density p = 1/4, using a 24 x 3 system with 
Nb = 18. The colored regions are dehneated using VMC 
data. Phase separation is determined with VMC methods by 
means of a Maxwell construction. Commensurate densities 
tend to be more stable for large K and so phase separation 
for this system generally consists of a region of zero density 
and a region at 1/3 density. DMRG points for the superfluid, 
DBL[3, 1], and phase separation are indicated by white circles, 
green squares, and gray diamonds, respectively. Our calcula- 
tions indicate that the DBL[3, 1] phase is never stabilized in 
the system with isotropic hopping (J±/J ~ 1). 



again we do not know whether the different signs of W2 
produce qualitatively distinct phases. 



Appendix B: Three- leg ladder at p 7^ 1/3 

In this paper, we have focused on the novel DBL[3,0] 
phase for the three-leg ladder at commensurate den- 
sity p = 1/3. We did explore, however, the non- 
commensurate densities on the three-leg ladder, and we 
shall summarize the results of those studies in this ap- 
pendix. 



1. DBL[3,1] at p = 1/4 < 1/3 



DMRG: 
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FIG. 15: (color online). Boson momentum distribution func- 
tion (top panel) and density-density structure factor (bottom 
panel) for a typical point in the DBL[3, 1] phase. The pa- 
rameters are J±/J — 0.75 and K/J = 2.5 on a system of 
size 24 X 3 with Nt = 18 (p = 1/4). DMRG data are joined 
with solid curves while the VMC data are joined with dashed 
curves. 



For densities smaller than the commensurate value, our 
data suggest that a stable DBL[3, 1] phase exists over a 
reasonable portion of the phase diagram (see Fig. 14). 
In particular, we studied p = 1/4 where the energy- 
optimized VMC state shows excellent qualitative agree- 
ment with the DMRG results. In the top panel of Fig.fTS] 
we show the comparison of the two data sets for the 
boson momentum distribution function while the bot- 
tom panel shows the density-density structure factor; the 
overall agreement is quite striking. The results are for a 
24 X 3 system with 18 bosons. The DBL[3, 1] state is con- 
structed as follows: the di partons partially fill all three 
bands with ten orbitals occupied for the ky — band and 
four orbitals occupied for each of the ky — ±27r/3 bands; 
meanwhile, the c?2 partons only fill the ky = band. The 
structure factors show enhanced features at several wave 



vectors as predicted by the corresponding gauge theory 
for this phase. Also, an analysis of the entanglement 
entropy with DMRG indeed reveals an effective central 
charge c ~ 3, as we expect from the gauge theory which 
leads to a theory ofc = 3 = 3 + l — 1 ID gapless modes. 

We do not give more attention to this phase since it is a 
rather straightforward extension of the DEL [2, 1] phase 
on the two-leg ladder that only exists for densities less 
than 1/A^, where N is the number of legs. The DBL[3, 1] 
on three legs, like the DBL[2, 1] on two legs, has exactly 
the same number of gapless modes as the number of legs 
and the c?2 determinant in the model wave function can 
be thought of as a Jordan- Wigner string, which is truly 
ID in nature. Furthermore, for = 4 at p < 1/4, we find 
that the system phase separates right out of the super- 
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fluid and no such DBL[4, 1] state is observed. Therefore, 
the DBL[3,1] is not extensible to 2D and is of limited 
interest. Nonetheless, its appearance is well supported 
by the underlying gauge theory giving further credence 
to our theoretical approach. 



2. Lack of DBL[n, m] phases for p > 1/3 

Here we summarize what we know about the situation 
for 1/3 < p < 1/2 on the three-leg ladder where we had 
hoped to find a DBL[3, 2] phase with one more gapless 
mode than the number of legs. In this putative phase, 
the di partons partially fill all three bands ky = 0, ±27r/3, 
whereas the d2 partons fully fill the ky = band and par- 
tially fill the ky = ±27r/3 bands. Such a phase would be 
extensible to 2D and constitute substantive progress in 
our DEL search. However, we performed an extensive 
study of the three-leg system at boson densities p = 1/2 
and 5/12, and no such phase was found. In fact, the VMC 
calculations do not even favor a DBL[3, 2] in any region 
of the phase diagram. Instead, with the VMC, we find 
a state in which the di partons do partially fill all three 
bands, while the c?2 partons fully fill the ky = band 
but partially fill only one of the remaining ky = ±27r/3 
bands; we denote this state DBL[3, 1]^ (with the dagger 
to distinguish it from the DBL[3, 1] discussed above in 
Appendix |B 1[ ) . Due to the population imbalance of the 
ky = ±27r/3 bands, this wave function breaks time re- 
versal symmetry and has a finit e b oson current in the 
y direction: (J^) ^ [see Eq. (18)] . This prompted 
us to scan the phase diagram with DMRG, consider- 
ing current-current correlations [see Eq. ( 19 )] in addition 



to the single-boson Green's function and density-density 
correlation function. 

We first focus on density p = 1/2, which is not com- 
mensurate on the three-leg ladder but, by particle-hole 
symmetry, is the largest possible boson density in our sys- 
tem. For simplicity, we also specialize our discussion here 
to J_l/ J = 1. While the VMC favors a DBL[3, 1]^ for in- 
termediate K (between a superfiuid at small K and phase 
separation at large K), we were able to rule out this phase 
with the DMRG. The state DBL[3, 1]"^ would have long- 
range, non-oscillatory y-y current-current correlations, 
i.e., Crir = (x,0),r' = (0,0)] ^ (jf (x, 0))(jf (0, 0)) ^ 
as a: — > oo, and a corresponding Bragg peak in C^'^(q) 
at q = 0. Interestingly, the DMRG measurement of 
Cj'^(q) for K/J > 2 shows clear Bragg peaks, not at 
q = 0, but at wave vectors on the corners of the Brillioun 
zone of our three-leg ladder: q = (tt, ±27r/3). There is 
also a weaker feature at q = (tt, 0). On the other hand, 
the boson momentum distribution function nf,(q) shows 
a distinct, superfiuid-like singularity at q = 0, as well as 
weaker features at q = (0, ±27r/3), (tt, ±27r/3). A set of 
DMRG measurements of nf,(q) and C^'^(q) for a char- 



acteristic point (Jl/J 
features is shown in Fig. 



1,K/J = 3) exhibiting these 
TE\ The density-density struc- 



ture factor (data not shown) is rather unremarkable: it 
shows (1) no sign of charge ordering and (2) \qx:\ behavior 
around g^; = at Qy = clearly indicating a compressible 
phase. Finally, for very large K, e.g., K/J = 10, we find 
that the three-leg half-filled system phase separates into 
p = 1/3 and p = 2/3 insulating regions. 

Although we have not attempted to understand 
this phase as an instability out of either DBL[3, 2] 
or DBL[3, 1]^, it does bear some resemblance to the 
phase, denoted "bond-chiral superfiuid" (BCSF), found 
recentljEH in a linear spin- wave treatment of the 2D J-K 
model [see Eq. this 2D phase is characterized by the 
coexistence of superfiuidity with static order in bond chi- 
rality (boson rung currents) at wave vector q = (tt, tt). 
To make this connection more concrete, we have per- 



formed a simple classical analysis of our model, Eq. (Ill, 



on the three-leg ladder, and we will now describe the 
results of this study. 

We first replace the boson operator with a classical c- 
number, i.e., — > y^exp(j0r), to obtain the classical 
energy: 



H — Hhop + Hk, 
"Hhop = - J^2pcos( 

r 

— Jl 2p cof. 

r 

Hk = K^2p^ cos( 



■'r+y h 



(Bl) 



(B2) 



Or+x + (Pr+x+y 



fr+y) 



(B3) 



Formally, this corresponds to relaxing the hard-core con- 
straint and computing the energy {H) = T-L with respect 
to a product of local coherent states with amplitudes 
Y^exp(i<^r)- Strictly speaking, such a state cannot ac- 
tually exist in our quasi-lD setting because it represents 
a superfiuid with true long-range order. However, as we 
describe next, on short scales a classical state with a par- 
ticular phase pattern 0r can qualitatively describe much 
of the DMRG data in Fig. [16) while the actual three- 
leg quantum phase would correspond to a slowly varying 
global rotation of the phase pattern on long scales. 



Numerical analysis of Eq. ( Bl ) on small periodic three- 



leg systems indicates that for substantial K the clas- 
sical ground state has an r-dependent phase pattern, 
4'r={x,y) — 4'y{x) (where y = 1, 2, 3), of the form 



-^1(2:) =0, ^2,3(2:) = ±(-i)= 



a . 



(B4) 



We depict this solution in the bottom panel of Fig. [16] 
by drawing at each site a vector (cos 0r, sin (/)r). This 
classical state has finite boson currents J^(r) = 
— 2/9sin (0r — 0r-i-/i) and thus long-range order in the 
current-current correlator Cg''*'(r, r'). The structure fac- 
tor C^'^(q) has Bragg peaks, i.e., Cjj'^{q)/Lx — 7 with 
7 a constant, at momenta q = (tt, ±27r/3), (tt, 0) with 
amplitudes 7 = |p^(sin a -I- sin 2q!)^, sin'* ^ sin^ a. 
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FIG. 16: (color online). Boson momentum distribution func- 
tion (top panel) and y-y current-current structure factor (mid- 
dle panel) for a typical point in the putative "bond-chiral su- 
perfluid" (BCSF) region at p = 1/2 on the three-leg ladder. 
The parameters are J±/J = 1 and K/J — 3 on two system 
sizes: Lx ~ 24, 12. As explained in the text, the features ap- 
pearing in these correlations can be qualitatively rationalized 
by a classical solution, br — ^/p exp {i<f)r) , with phases given 



by Eq. (B4 1. In the bottom panel, we display this solution by 
showing a vector (cos^r, sin(/!>r) at each site. Since true long- 
range order of the boson phase is prohibited for our three-leg 
ladder, this is a potentially appropriate description only on 
short length scales. Note, however, that the current-current 
correlations can still be truly long ranged. 



respectively. To get a very rough idea of the num- 
bers involved, Eq. (Bl) is minimized at p = 1/2 and 
Jj_/J = 1, K/J = 3 by Eq. ([B4| with a = 0.73, which 
gives Bragg peaks at q = (tt, ±27r/3), (tt, 0) with respec- 
tive amplitudes 7 = 0.92,0.04. This is in quahtative 
agreement with the apparent large Bragg peak in the 
DMRG data at (7r,±27r/3) but only a very weak fea- 
ture at (7r,0) (see Fig. 16 middle panel). (Note that 



the ratio K/ J is even somewhat arbitrary in this anal- 
ysis since it is affected by the normalization of the bo- 
son field, which we take as 6j6r = p, cf. Rcf. 33 which 
formally takes &j6r = 1/4 = p/2.) The classical solu- 
tion also has a Bragg peak in the x-x current-current 
structure factor C^'^(q) at q = (tt, ±27r/3) with ampli- 
tude 7 = sin^ 2a. We have measured x-x current- 
current correlations with the DMRG (data not shown), 
and although they appear weaker than the correspond- 
ing y-y correlations, we still cannot rule out a Bragg 
peak at (tt, ±27r/3). Finally, one can compute the bo- 
son momentum distribution within this naive framework. 
This calculation predicts Bose condensation at wave vec- 
tors q = (0, 0), (0, ±27r/3), (tt, ±27r/3), which are inter- 
estingly the same wave vectors at which we see fea- 
tures in the DMRG. Specifically, at these wave vec- 
tors, rtf,(q)/La; = |/9(1 -1-2 cosa)^, |yOsin'' ^, psin^ a , re- 
spectively. To again get a feel for the numbers, these 
three amplitudes evaluated at p = 1/2 and a = 0.73 
give 1.03,0.01,0.22. Again, this simple classical anal- 
ysis agrees fairly well qualitatively with the DMRG, 
which indeed indicates possible quasi-condensation at 
momenta q = (0, 0), (0, ±27r/3), (tt, ±27r/3) (see Fig. [I6 
top panel). 

The fact that such an exceedingly crude classical anal- 
ysis qualitatively reproduces much of the DMRG data is 
quite remarkable. However, we caution that the above 
picture is likely not complete and could still be qualita- 
tively wrong. Indeed, more work can certainly be done 
to characterize this unusual quantum phase. For exam- 
ple, we have been unable to extract a meaningful central 
charge with the DMRG data, which points to c ~ 2-3 for 
a 24x3 system, while the quasi-lD BCSF phase described 
above should lead to a Luttinger liquid with c — 1 gapless 
mode. It may be that the DMRG entanglement entropy 
is anomalously large due to the large ground state de- 
generacy in the putative BCSF phase, and c = 1 scaling 
only emerges at much longer length scales. 

Upon doping the half-filled three-leg system to boson 
density p = 5/12, the putative BCSF phase discussed 
above appears to persist, at least for sizable J± and in- 
termediate K, e.g., J±/J = 1 and 2 < K/J < 3.5. For 
smaller values of J±, e.g., = J±/J ~ 0.1, we rather re- 
markably find DMRG evidence for the DBL[3, 1]^ state 
in which the di partons equally fill the three bands. In 
this case, there is a conserved number of bosons in each 
chain and thus zero rung current in the projected wave 
function: (J^) = 0. Loosely speaking, the small value 
of J± has a tendency to suppress rung currents and 
seems to stabilize the decoupled chains limit of so-called 
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DBL[3,1]''', though this is the only incarnation of this 
phase that we have been able to realize in the DMRG. 
Also, DBL[3, 1]^ does not appear to exist at p = 1/2, even 
at small Jl, which is at least somewhat surprising. Be- 
cause the DBL[3, 1]^ is obviously very far removed from 
2D, we do not pursue understanding it any further. 

In conclusion, the identified putative three-leg BCSF 
phase (see Fig. 16 1 is itself rather interesting and wor- 



thy of further investigation. We stress that we found no 
evidence for such a phase on the four-leg ladder at any 
density; e.g., we found no evidence of long-range order in 
current-current correlations. In that case, our DBL[4, 2] 
is stabilized for densities slightly below p = 1/2 (e.g., 
p — 5/12), while ai p — 1/2, the conventional q = 
superfluid was the only compressible phase identified. 
Thus, exactly how all of these p > 1/3 th ree-l eg results 
and our p > 1/4 four-leg results (see Sec. Ill) extrapo- 



late to 2D, if at all, is not presently clear. In Sec. Ill 



we presented strong evidence in support of the existence 
and stability of a very 2D-like DBL phase on four legs, 
namely, the so-called DBL[4, 2]; however, as we have dis- 
cussed in this appendix, an analogous phase does not 
appear to exist for the three-leg system. Generally speak- 
ing, we do believe odd-leg systems are somewhat patho- 
logical when approaching 2D in this context because they 
are not bipartite, thus putting the half-filled system on 
qualitatively different footing than what we would ex- 
pect in 2D. Still, to further investigate this seemingly 
strong even-odd effect, in addition to going to more legs 
(which is becoming exceedingly difficult at this point) , it 
may be insightful to modify the transverse (y direction) 
boundary condition for the boson model [Eq. (|TT|)]. In 
this work, we have only considered periodic transverse 
boundary conditions in both the three- and four-leg sys- 
tems. However, the structure of the DBL[4, 2] identi- 
fied on four legs suggests that an anti-periodic boundary 
condition may be desirable to stabilize the sought-after 
DBL[3, 2] for p > 1/3 on three legs. Namely, because the 
structure of the stable four-leg DBL[4, 2] relies on anti- 
periodic transverse boundary conditions for both the di 
and d2 partons, a more natural three-leg DBL[3,2] con- 
figuration than that discussed at the beginning of this 
appendix might consist of a periodic boundary condition 
for di and an anti-periodic boundary condition for d2. 
This configuration is only possible with an anti-periodic 
boundary condition for the boson system. We leave such 
further exploration for future work. 

These effects illustrate the degree of difficulty encoun- 
tered when trying to extrapolate results from quasi-lD 
ladder studies to make definitive statements about 2D 
physics, especially when trying to understand a phase as 
complicated as the DBL in a nontrivial model such as 
the J-K model. However, we believe our four-leg results 
of Sec. |IIIB| are more representative of the 2D system 
than the three-leg results summarized above in this ap- 
pendix: firstly, in contrast to the three-leg ladder, the 
four-leg ladder is a bipartite lattice (also, four is greater 
than three), and secondly, fully periodic boundary con- 



ditions are the most natural choice for approaching full 
2D. 



Appendix C: Two- leg ladder at p = 1/2 

We now address the situation on the two-leg ladder 
at half-filling (p = 1/2). In this system, it is in prin- 
ciple possible for a gapless Mott insulating phase, i.e., 
a DBL[2,0] in our naming convention, to exist for large 
K . A nontrivial solvable limit is when J = 0, in which 
the number of bosons in each rung is conserved and is 
exactly one for p = 1/2. In this regime, we can map the 
model exactly to the spin- 1/2 XY chain with an in-plane 
magnetic field (identifying a boson on the upper chain as 
spin up and a boson on the lower chain as spin down): 



H 



XY 



i+l 



H.c. 



(Cl) 



where a^''^'^ are the Pauli matrices for the pseudospin 
on rung i. When Ji — 0, the XY model is solved by 
a Jordan- Wigner transformation and has one ID gapless 
mode corresponding to the free fcrmion. Jj^ is relevant at 
this XY fixed point and immediately causes the forma- 
tion of a canted magnet (uniform moment along x and 
staggered moment along y) before entering the trivial 
pseudospin paramagnet (staggered moment along y goes 
to zero) at large Jj^.^-'^ In the two-leg ladder boson lan- 
guage, this corresponds to a staggered rung current state 
and the featureless rung Mott insulator, respectively; the 
transition between these two gapped phases is continu- 
ous and is of the Ising universality class. The phase dia- 
gram at J = can be logically extended into the entire 
K/ J ~ J_l/ J plane without the introduction of an exotic 
GMI phase. Indeed, we have verified this with DMRG 
and found no evidence for a GMI phase at nonzero J_l . In 
addition to the staggered rung current phase mentioned 
above, a (tt, tt) CDW appears at intermediate K/ J and 
small Ja_/J. Remarkably, a bosonization analysis of the 
putative two-leg GMI phase predicts an allowed insta- 
bility toward either (tt, tt) static order in the density or 
TT static order in the rung currents, both of which we 

actually observe in the DMRG. 

Finally, we note that two very recent workP^l^ have 
considered in detail hard-core bosons hopping on the two- 
leg ladder without ring exchange, i.e., our J-J±-K model 
with K = 0. In these works, it was shown that at half- 
filling the system enters a rung Mott insulating phase for 
any finite Jl, although the charge gap grows exponen- 
tially slowly with Jj_, i.e., as exp(— a J/ Jj^), and is ex- 
tremely difficult to deal with numerically due to a large 
value of the constant factor a in the exponential.'^ In 
our DMRG study of the two-leg half-filled system, we ob- 
served what looked like a superfluid on finite-size systems 
for small K; however, we did not perform an extensive 
finite-size analysis of this system as in Ref . [56l It is there- 
fore possible that the apparent superfluid is a rung Mott 
phase on long length scales, and there are only direct 
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transitions between the rung Mott phase and the stag- 
gered rung current and (7r,7r) CDW phases mentioned 
above. 

To conclude, although we were able to successfully 
identify a GMI phase on the three-leg ladder (see Sec. IV 



and Ref . [TO]) , analogous phases do not exist on the two- 
and four-leg systems, at least with our model, boundary 
conditions, etc. This again suggests a rather strong even- 
odd effect similar to what we found above for conducting 
DBL phases at incommensurate densities. 
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